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ABSTRACT. We consider a family of discrete coagulation-fragmentation equa- 
tions closely related to the one-dimensional forest-fire model of statistical me- 
chanics: each pair of particles with masses i,j 6 N merge together at rate 2 to 
produce a single particle with mass i + j, and each particle with mass i breaks 
into i particles with mass 1 at rate (i — l)/n. The (large) parameter n controls 
the rate of ignition and there is also an acceleration factor (depending on the 
total number of particles) in front of the coagulation term. We prove that for 
each n £ N, such a model has a unique equilibrium state and study in details 
the asymptotics of this equilibrium as n — > oo: (I) the distribution of the mass 
of a typical particle goes to the law of the number of leaves of a critical binary 
Galton- Watson tree, (II) the distribution of the mass of a typical size-biased 
particle converges, after rescaling, to a limit profile, which we write explicitly 
in terms of the zeroes of the Airy function and its derivative. We also indicate 
how to simulate perfectly a typical particle and a size-biased typical particle, 
which allows us to give some probabilistic interpretations of the above results 
in terms of pruned Galton- Watson trees and pruned continuum random trees. 
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1. Introduction 

1.1. The forest-fire model. The forest-fire model of statistical mechanics has 
been introduced by Henley [25] and Drosscl-Schwabl [16] in the context of self- 
organized criticality. From the rigorous point of view, the one-dimensional forest- 
fire model has been studied by Van den Berg-Jarai [8], Brouwer-Pcnnanen [12j and 
by the authors [TDJ [TT] . We refer to the introduction of [TT] for many details. 

Let us now describe the one-dimensional forest-fire model: on each site of Z, 
seeds fall at rate 1 and matches fall at rate 1/n, for some (large) neN. Each time 
a seed falls on a vacant site, this site immediately becomes occupied (by a tree). 
Each time a match falls on an occupied site, it immediately burns the corresponding 
occupied connected component. 

From the point of view of self-organized criticality, it is interesting to study what 
happens when n increases to infinity. Then matches are very rare, but tree clusters 
are huge before they burn. In |10| . we have established that after normalization, 
the forest-fire process converges, as n — > oo, to a scaling limit. 

1.2. A related mean-field model. We now introduce a mean-field model for- 
mally related to the forest-fire process, see [9] Section 6] for a similar study when 
n = 1. Assume that each edge of Z has mass 1. Say that two adjacent edges (i — 1, i) 
and (i, i + 1) are glued if the site i is occupied. Then adjacent clusters coalesce at 
rate 1 and each cluster with mass k (containing k edges and k — 1 sites) breaks up 
into k clusters with mass 1 at rate (k — l)/n. Denote by du(t) the concentration 
(number per unit of length) of clusters with mass k > 1 at time t > 0. Then 
the total mass should satisfy J2k>i ^' C L- W = 1 f° r a ^ t > 0. Neglecting correla- 
tions (which is far from being justified), the family (c£ (i))t>o,fc>i would satisfy the 
following system of differential equations called (CF n ): 

(1) ±$(t) = 2c?(t) + - ]T k(k - l)cl(t), 
at n f-^ 

k>2 

(2) £%(t) = (2 + (k - l)/n)ct{t) + \ J2 c t(t)cU(t) (k > 2). 

The first term on the RHS of ([1} expresses a cluster with mass 1 disappears at rate 
2 (because it glues with each of its two neighbors at rate 1); the second term on 
the RHS of ([1]) says that k clusters with mass 1 appear each time a cluster with 
mass k takes fire, which occurs at rate (k — l)/n. The first term on the RHS of 



A MEAN FIELD FOREST-FIRE MODEL 



3 



([2]) explains that a cluster with mass k disappears at rate 2 + (k — l)/n (because 
it glues with each of its two neighbors at rate 1 and takes fire at rate (k — l)/n). 
Finally, the second term on the RHS of @ says that when a seed falls between two 
clusters with masses i and k — i, a cluster with mass k appears. The number per 
unit of length of pairs of neighbor clusters with masses i and k — i is nothing but 
c ™(^) c fe-iW/ X)i>i C FW- Here we implicitly use an independence argument which 
is not valid for the true forest-fire model. 

The system (d])-© can almost be seen as a special coagulation-fragmentation 
equation, see e.g. Aizcnman-Bak [2] and Carr |13j . where particles with masses i 
and j coalesce at constant rate K(i,j) = 2 and where particles with mass i > 2 
break up into i particles with mass 1 at rate F(i; 1, . . . , 1) = (i — l)/n. However 
there is the acceleration factor 1/Y^i>i C "W m front of the coagulation term. 

1.3. On the link between the two models. The link between (CF n ) and the 
forest-fire model is only formal. Observe however that if there are only seeds or 
only matches, then the link is rigorous. 

(i) Assume that all sites of Z are initially vacant and that seeds fall on each site 
of Z at rate 1, so that each site is occupied at time t with probability 1 — e~*. Call 
Pk(t) the probability that the edge (0, 1) belongs to a cluster with mass k at time t. 
A simple computation shows that pk(t) = k(l — e~ t ) k ~~ 1 e~ 2t . By space stationarity, 
the concentration Ck(t) of particles with mass k at time t satisfies Ck(t) = Pk(t)/k = 
(1 — e~ t ) k ~ 1 e~ 2t . Then one easily checks that the family (ck(t))k>i,t>o satisfies 
CE])-© with no fragmentation term (i.e. n = oo). 

(ii) The fragmentation term is linear and generates a priori no correlation. As- 
sume that we have initially some given concentrations (cfc(0))fc>i, so that the edge 
(0, 1) belongs to a cluster with mass k with probability kck(0) and that only matches 
fall, at rate 1/n on each site. Then the probability pk(t) that the edge (0, 1) belongs 
to a cluster with mass k at time t is simply given by Pk(t) = fccfc(0)e - ( fc-1 - ) */™ if 
k > 2 (here e-^- 1 )*/™ i s the probability that no match has fallen on our cluster be- 
fore time t) andpi(t) = ci(0) + £ fc>2 fcc fe (0)[l- e"^" 1)t/n ] . Writing c k (t) =p k (t)/k 
as previously, we see that the family (ck(t))k>i,t>o satisfies the fragmentation equa- 
tions |ci(t) = n~ x Y, k > 2 k{k- l)cfc(t) and, for k > 2, f t c k (t) = - n -\k - l)c k (t). 

When one takes into account both coalescence and fragmentation, the rigorous 
link between the two models breaks down: in the true forest-fire model, fragmen- 
tation (fires) produces small clusters which are close to each other, so that a small 
cluster has more chance to have small clusters as coalescence partners. However, 
we have seen numerically in Section 6] that at equilibrium, in the special case 
where n = 1, the two models are very close to each other. 

1.4. Motivation. Initially, the motivation of the present study was to decide if (p}- 
(J2j) is a good approximation of the true forest-fire model, at least from a qualitative 
point of view: do we have the same scales and same features (as n — > oo)? We 
will see that this is not really the case. However, we believe that (HJ-© is a very 
interesting model, at least theoretically, since 

• many explicit computations are possible for CO)-©, 

• we observe self-organized criticality, 

• we show that two interesting points of view (size-biased and non size-biased 
particles' mass distribution) lead to quite different conclusions, 
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• we have some probabilistic interpretations of our results. 

1.5. Summary of the main results of the paper. We will show in this paper 
the existence of a unique equilibrium state (c^)k>i with total mass J2k>i ^ c fe = 1 
for (CF n ), for each n > 1 fixed and we study the asymptotics of rare fires moo. 

(I) We show that the particles' mass distribution (P/J)fc>i, defined by p£ = 
c k/ J2i>i c Y' S oes weakly to the law {pk)k>i of the number of leaves of a critical 
binary Galton- Watson tree, which is explicit and satisfies pk ~ (2 y /nk 3 ^ 2 )~ 1 as 
k — > oo; 

(II) We prove that the size-biased particles' mass distribution (fcc£)fc>i goes 
weakly to, after normalization of the masses by n~ 2 / 3 , to a continuous limit profile 
(xc(a;))a; e (o,oo) with total mass 1, for which we have two explicit expressions: the 
first one involves the zeroes of the Airy function and its derivative, while the Laplace 
transform of the Brownian excursion's area appears in the second one. Furthermore, 
we check that c(x) ~ K\x~ z l 2 as x — > and c(x) ~ K2e~ K3X as x — > oo, for some 
positive explicit constants Ki,K2,K3- 

(III) We explain how to simulate perfectly, for n > 1 fixed, a random variable X n 
with law (pk)k>i an d a random variable Y n with law (kc%.)k>i, using some pruned 
Galton- Watson trees. 

(IV) We deduce from (III) an easy probabilistic interpretation of point (I): 
roughly, (p%)k>i can be viewed as the law of the number of leaves of a pruned 
critical Galton- Watson tree, and we only have to check that almost no pruning 
occurs when n is very large. 

(V) We also derive a probabilistic interpretation of (II), which seems quite nat- 
ural, since the Brownian excursion's area appears in the limit. This part is quite 
complicated. The main idea is that that the pruned Galton- Watson tree of which 
the number of leaves is (fcc£)fc>i-distributed has a scaling limit, which is nothing 
but a pruned continuum random tree (precisely, a pruned self-similar CRT, see 
Aldous [3]). We do not prove rigorously that this scaling limit arises. However, we 
show that in some precise sense, the number of leaves of the pruned CRT under 
consideration is indeed a;c(x)(ia;-distributed, where c is the limit profile found in 
point (II). 

(VI) Finally, this pruned CRT (heuristically) leads to a noticeable diffusion pro- 
cess (rigorously) enjoying the strange property that its drift coefficient equals the 
Laplace exponent of its inverse local time at 0. 

Let us discuss briefly points (I) and (II). The particles' mass distribution is, 
roughly, the law of the mass of a particle chosen uniformly at random. The size- 
biased particles' mass distribution is, roughly, the law of the mass of the particle 
containing a given atom, this atom being chosen uniformly at random (think that 
a particle with mass k is composed of k atoms). 

Point (I) says that if one picks a particle at random, then its mass X n is finite 
(uniformly in n) and goes in law, as n — > oo, to a critical probability distribution 
(with infinite expectation). Point (II) says that if one picks an atom at random, 
then the mass Y n of the particle including this atom is of order n 2 / 3, and n~ 2 l 3 Y n 
goes in law to an explicit probability distribution with moments of all orders. 
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1.6. Comments. Let us now comment on these results. 

(a) Self-organized criticality, see Bak-Tang-Wicsenfield [7] and Henley [35], is a 
popular concept in physics. The main idea is the following: in statistical mechanics, 
there arc often some critical parameters, for which special features occur. Consider 
e.g. the case of percolation in Z 2 , see Grimmett [21]: for p £ (0, 1) fixed, open each 
edge independently with probability p. There is a critical parameter p c G (0, 1) 
such that for p < p c , there is a.s. no infinite open path, while for p > p c , there is 
a.s. one infinite open path. If p < p c , all the open paths are small (the cluster-size 
distribution has an exponential decay). If p > p c , there is only one infinite open 
path, which is huge, and all the finite open paths are small. But if p = p c , there 
are some large finite open paths (with a heavy tail distribution). And it seems that 
such phenomena, reminiscent from criticality, sometimes occur in nature, where 
nobody is here to finely tune the parameters. Hence one looks for models in which 
criticality occurs naturally. We refer to the introduction of |11) for many details. 

(b) Thus we observe here self-organized criticality for the particles' mass distri- 
bution, since the limit distribution (pk)k>i has a heavy tail and is indeed related 
to a critical binary Galton- Watson process. Observe that point (I) is quite strange 
at first glance. Indeed, when n — > oo, the time-dependent equations CQ)-© teno - to 
some coagulation equations without fragmentation, for which there is no equilibrium 
(because all the particles' masses tend to infinity). However, the equilibrium state 
for (HJ-© tends to some non-trivial equilibrium state as n — > oo. This is quite sur- 
prising: the limit (as n — > oo) of the equilibrium is not the equilibrium of the limit. 
Observe that this is indeed self- organized criticality: the critical Galton- Watson 
tree appears automatically in the limit n — > oo. A possible heuristic argument is 
that in some sense, in the limit n — > oo, only infinite clusters are destroyed. We 
thus let clusters grow as much as they want, but we do not let them become infinite. 
Thus the system reaches by itself a critical state. 

(c) For the size-biased particles' mass distribution, we observe no self-organized 
criticality, since the limit profile has an exponential decay. The two points of view 
(size-biased or not) seem interesting and our results show that they really enjoy 
different features. Let us insist on the fact that the particles' mass distribution 
converges without rescaling, while the size-biased particles' mass distribution con- 
verges after rescaling. This is due to the fact that there are many small particles 
and very few very large particles, so that when one picks a particle at random, we 
get a rather small particle, while when one picks an atom at random, it belongs to 
a rather large particle. Mathematically, since the limiting particles' mass distribu- 
tion has no expectation, it is no more possible to write properly the corresponding 
size-biased distribution: a normalization is necessary. 

(d) Let us mention the paper of Rath-Toth [32] , who consider a forest-fire model 
on the complete graph. In a suitable regime, they obtain the same critical distri- 
bution (pk)k>i as w e do (see [32] Formula (14)]). But in their case, this is the limit 
of the size-biased particles' mass distribution. Their model rather corresponds to 
the case of a multiplicative coagulation kernel (clusters of masses fc and I coalesce 
at rate kl). Hence we exhibit, in some sense, a link between the Smoluchowski 
equation with constant kernel and the Smoluchowski equation with multiplicative 
kernel. See Remark [3] for a precise statement. In the same spirit, recall the link 
found by Deaconu-Tanre |15j between multiplicative and additive coalescence. 
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(e) For the true forest-fire process, we proved in |10| the presence of macroscopic 
clusters, with masses of order n/logn and of microscopic clusters, with masses 
of order n z , for all values of z G [0, 1). These microscopic clusters interact with 
macroscopic clusters in that they limit the impact of fires. Hence, while the scales 
are really different here, we observe something similar: in some sense, the particles' 
mass distribution describes microscopic clusters, while the size-biased the particles' 
mass distribution describes macroscopic clusters. 

(f) The trend to equilibrium for coagulation-fragmentation equations has been 
much studied, see e.g. Aizenman-Bak [2] and Carr [13], under a reversibility condi- 
tion, often called detailed balance condition. Such a reversibility assumption cannot 
hold here, because particles merge by pairs and break into an arbitrary large num- 
ber of smaller clusters. Without reversibility, much less is known: a special case has 
been studied by Dubowski-Stcwart |17j and a general result has been obtained in 
[21] under a smallncss condition saying that fragmentation is much stronger than 
coalescence. None of the above results may apply to the present model (at least 
for n slightly large). For the specific model under study we are able to prove the 
uniqueness of the equilibrium, but not the trend to equilibrium. 

(g) Finally, let us mention that our probabilistic interpretation of (II) uses ex- 
tensively Aldous's theory on continuum random trees [3] EJ [5], which has been 
developed by Duquesne-Le Gall [TH [20] . Some pruning procedures of such CRTs 
have been introduced by Abraham-Dclmas-Voisin [T] . Roughly, the prune the CRT 
by choosing cut-points uniformly in the tree (uniformly according to the Lebcsgue 
measure on the tree). Our pruning procedure is very different: we choose leaves 
(according to the measure on leaves), these leaves send some cut-points on the 
branches (joining them to the root). Then we prune according to these cut-points, 
in a suitable order. The cut-points sent by leaves belonging to subtrees previously 
pruned are deactivated. 

1.7. Outline of the paper. The next section is devoted to the precise statement 
of points (I) and (II). Section [3] contains the proofs of points (I) and (II), which are 
purely analytic. Using discrete pruned random trees, we indicate how to simulate 
perfectly a typical particle and a size-biased typical particle in Section 2] Finally, 
we study in Section [5] a possible scaling limit of the discrete tree representing a 
typical size-biased particle. 



Definition 1. Let n 6 N. A sequence c" = (c^)k>i of nonnegative real numbers 
is said solve {E n ) if it is an equilibrium state for (CF n ) with total mass 1: 



2. Precise statements of the results 



For a [0, oo)-valued sequence u = (uk)k>i and for i > 0, we put 




(3) m x {c n ) =1, 
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Observe that it then automatically holds that 

( 5) 2c; = i£ M4 _ 1K = _*!) 



n 

k>2 



To check this last claim, multiply (Q| by k, sum for k > 2, use © and that 

£ fc > 2 fcEti 1 = E fc > 2 E-=i(i + * " *K<£_< = 2E M >! fccjcj* = 2m„(c»). 

One finds 2(1 — c") + [m2(c n ) — l]/n = 2, from which ([S]) readily follows. 

To state our main results, we need some background on the Airy function Ai. 
We refer to Janson [2j)J P 94] and the references therein. Recall that for x £ R, 
Ai(x) = tt" 1 f o °° cos(i 3 /3 + xt)dt is the unique solution, up to normalization, to 
the differential equation Ai"(x) = xA.i{x) that is bounded for x > 0. It extends 
to an entire function. All the zeroes of the Airy function and its derivative lie on 
negative real axis. Let us denote by a[ < the largest negative zero of Ai' and by 
• • • < 03 < 02 < ai < the ordered zeroes of Ai. We know that |ai| ~ 2.338 and 
\a' x \ ~ 1.019, see Finch [23]. We also know that \a,j\ ~ (37r?'/2) 2 / 3 as j -> 00, see 
Janson [H p 94]. 

Theorem 2. (i) For each n € N, (E n ) has a unique solution c n = (c2)k>i- 
(ii) As n — > 00, there hold 

1 n 2 / 3 

In some sense, mo(c") stands for the total concentration and 7712(0") stands for 
the mean mass of (size-biased) clusters. For any fixed / > 0, we also have shown 
that 7Tii + i(c") ~ Min 2l / n for some positive constant Mi, see Lemma [9] below. 

Theorem 3. For each n > 1, consider the unique solution {d£)k>l to (E n ) 
and the corresponding particles' mass probability distribution (p^)k>i defined by 
Pk = c fc/ TO o(c n ). There holds 

lim o ^2\pt- Pk \ = 0, where Pk :=JLf 2 *-* 

k>l v 

The sequence (pk)k>i is the unique nonnegative solution to 

^ k-l 

(6) = L Pk = -^^PiPk-i (k>2). 

k>l i=l 
There holds, as k — > 00, 

(7) Pk 



2v^t 3 / 2 ' 

Formally, divide ((4]) by mo(c") and make n tend to infinity: one gets 2 P k = 
Ej=i PiPk-% f° r au A: > 2. What is much more difficult (and quite surpris- 
ing) is to establish that no mass is lost at the limit. Observe that ^ rewrites 

Ei>iPi-Pfc = (1/2) Ei=i PiPk-i (for all k > 2), which corresponds to an equilib- 
rium for a coagulation equation with constant kernel. This is quite strange, since 
coagulation is a monotonic process, for which no equilibrium should exist. The 
point is that in some sense, infinite particles are broken into particles with mass 1, 
in such a way that J2k>iPk = 1- Finally, we mention that (pk)k>i is the law of the 
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number of leaves of a critical binary Galton- Watson tree, which will be interpreted 
in Section |U 

Theorem 4. For each n > 1, consider the unique solution (c2)k>i to (E n ) and 
the corresponding size-biased particles' mass probability distribution (kc2)k>i- For 
any (f> € C([0,oo)) with at most polynomial growth, there holds 



lim V (j)(n- 2/3 k)kc% = / 
fe>l • /u 



4>{x)xc(x)dx, 



where the profile c : (0, oo) i— > (0, oo) is defined, for x > 0, by 

oo 

c(x) = |ai| _1 exp(|ai|a;)^exp(-|aj|a:) . 

j'=i 

The profile c is of class C°° on (0, oo), has total mass / °° xc(x)dx = 1 and 

c O) X ~° i- x - i 3/2 an d c(x) x ~°° |ai| _1 exp((|a' 1 | - \ai\)x) . 
2y'Tr\a 1 \x- i ' z 

For any <f> € C 1 ([0,oo)) such that <fi and <fi' have at most polynomial growth, 

(8) 2\a' 1 \ / x[<j){x + y) - (j>(x)]c(x)c(y)dydx = / x 2 [4>(x) - (f>(0)]c(x)dx. 
Jo Jo Jo 

Denote by B ex is the integral of the normalized Brownian excursion, see Revuz- Yor 
pi Chapter XII] . For all x > 0, 



cx P (| ai |x) 
1 ' 2V¥|a' 1 |x3/2^ 



Finally, for all q 6 (ai — a' l7 oo) (recall that a\ — a[ < 0), 
(9) %) := f°(l - e-'-JcC^di - ~ Ai ' (<Z + a ' l] 



|ai|Ai(g + oi) 

Since the mean mass of a typical (size-biased) particle is of order n 2 / 3 by Theorem 
H}(ii), h is natural to rescale the particles' masses by a factor rt~ 2 / 3 . Here we state 
that under this scale, there is indeed a limit profile and we give some information 
about this profile. The link with the Brownian excursion's area will interpreted in 
Section [5j 

3. Analytic proofs 
For each n € N, we introduce the sequence (ct^)k>i, defined recursively by 

k-i 

(10) a? = l, (2 + (fc-l)/n) = (k > 2). 

i=i 

We also introduce its generating function /„, defined for q > by 

(11) /„(<?) = £>£g fe . 

k>l 

Obviously, /„ is increasing on [0, oo) and takes its values in [0, oo) U {oo}. The 
unique solution to (E n ) can be expressed in terms of this sequence. 
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Lemma 5. Fix neN. Assume that there exists a (necessarily unique) q n > 
such that f n {q n ) = 1. Assume furthermore that f' n (q n ) < oo. Then there is a 
unique solution to {E n ) and it is given by 

(12) $ = c$fi- 1 /f n {q n ) (fc>l). 

Furthermore, there holds mo(c") = l/(q n f^(q n )). 

Proof. We break the proof into 3 steps. 

Step 1. A simple computation shows that for any fixed r > 0, x > 0, the sequence 
defined recursively by 

^ fc-i 

yi = x, (2 + (fc - l)/n) j/fc = - y] yiVk-i (k>2) 

i=l 

is given by yk = a^x(x/r) 1 . 

Step 2. Consider a solution (c^)k>i to (E n ). Then due to ((3]) and Step 1 (write 
r = too(c") and x = c"), cJJ = a 7 k l c"(ci /m (c n )) k ^ 1 . We deduce that 

1 = £ c * = E «feKMo(c")) fe = / n (^/mo(0). 

m ° lc j fe>i fc>i 

Consequently, c™/m (c") = g„, whence c£ = a^c"g^ _1 . Next we know that 
mi(c") = 1, so that fn{q n ) = J2k>i k° t k ( ln~ 1 = mi(c n )/cj* = 1/c". Consequently, 
c i = l//4(?n) and thus c'fc = oQq*' 1 / f' n (q n ) as desired. 

5iep 5. Let us finally check that c™ as defined by (fT2|) is indeed solution to (E n ) 
and that it satisfies mo(c n ) — l/(q n f' n (q n )). First, 

/ n\ 1 n fe-i fn\1n) 1 



]g ^ Mqn) qnf' n {qny 

Next, © holds, since 

777, (V'M - 1 V hr, n n k - 1 - ^ qn " > - I 

mi(c ~ _L 

Rewriting c™ = a^a;(a;/r) fe - :L with x = 1/ f' n {q n ) and r = l/{q n f n {q n )), Step 1 
implies that for k > 2, 

fe-l fe-i 
(2 + (fc - l)/„) eg = q n f' n {q n ) £ 0?CU = —r- £ 

whence Q. □ 

To go on, we need some background on Bcssel functions of the first kind. Recall 
that for fc > and z € C, 



(13) = |r £ 



(-l)'z 



2 fe ^ 4^!(/ + fc)!' 
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We have the following recurrence relations, see [6j Section 4.6]: for k > 0, x G K, 



(14) 
(15) 
(16) 



Jfc+ 2 (x) 



2{k + l) 



J k +i(x) - Jfc(x), 



.k+I 



Jfc(x). 



It is known, see [51 Section 4.14], that all the zeroes of Jk are real. For all k > 1, 
we denote by jk the first positive zero of J k and by j' k the first zero of J' k . The 
sequence {jk)k>o is increasing, see Section 4.14]. Furthermore, < j' k < j k for 
all k > 1 , see [35J page 3] . We will also use that there exists a constant C > such 
that for all k > 1, see [22] pages 2-3], 

(17) k + 2- 1 / 3 |a' 1 |fc 1 /3 < j' k <k + 2- 1 /3| fl ' i | fc l/3 + Ck ~l/3 } 

where a! x is, as previously defined, the largest negative zero of Ai'. 

Lemma 6. Let n 6 N be fixed. The radius of convergence of the entire series 
fn is r n = (j2„_ i/(2n)) J2. For all x £ [0,r n ), there holds 



(18) 



fn(x) 



'2x 



J 2n (2nV2x) 



J 2 „-i(2nV2x) 

There exists a unique q n € (0,r„) such that f n (q n ) = 1< There holds f' n {q n ) 
(n/q n )(2q n — 1 + l/?i) € (0, oo). Finally, as n — > oo, 

1 r 



1 + Kln" 2 / 3 +0(n- 1 ). 



Proof. We fix n G N and divide the proof into five steps. 

Step 1. Put s„ = (j2n-i/(2?i)) 2 /2. By definition of jin-\ and since J2n-i is odd 
on R and has no complex zeroes, J 2n _i (2ny/2z) does not vanish on {0 < |z| < s„} C 
C. We thus may define g n (z) := s/2zJ2 n (2ns/2z) / J2 n -i(2ny/2z) on {0 < |z| < s„}. 
Using (|13[) . one easily checks that, as z —> 0, g n (z) ~ z, so that finally, g n (z) is 
holomorphic on the disc {\z\ < s n }. Write g n (z) = J2 k > f3 k z k . Since g n {z) ~ z 
near 0, we deduce that (3^ = and /3" = 1. 

Siep ,2. We now show that for all x G [0, s„), there holds 
(19) xg' n {x)/n + (2 - l/n)<? n (x) = g 2 (x) + 2x. 

Write g n {x) = h n (2n\/2x) / \2n) , where h n (y) = yhniv)/ '<^2n-i(j/)- Observing that 
MlO = (y 2 "^2n(2/))/(y 2 ™" 1 J 2n _i(y)) and using JTBJ and then (JTSJ, 

2/ 2 "J 2n -i(2/) 2/ 2 ™J2, i (y)2; 2 "- 1 J 2n _ 2 ( 2 /) 



y^J^-M (y 2n - 1 J2n-i(y)) 2 

J2n-2(y) 



=y- K{y) 
=y - K{y) 



J2n-l(y) 

-J 2n (y) | 2(2n-l) 
J2n-i(y) y 
_ 2(2n _l)M/) 
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But g' n (x) = h' n (2nv2x)/v2x, whence 

X JM. = * \ 2 nV^ + h ^ 2n ^3 - 2(2n - 1)^^ 
n nV2x [ 2nV2x 2nV2x 

=2x + gl{x)-(2-l/n)g n {x). 

Step 3. Let us check that the sequence (/3^ )fc>i satisfies (|10[) . We already know 
that P? = l. Using (|19|. 

fc-i 

£ WO** + ^(2 - l/n)/3^ fe = W-< + 

fe>l *;>1 k>2 t=l 

Thus for all fc > 2, (fc/n)j9£+(2-l/n)j9£ = ^=1 as desired. Consequently, 

(Pk)k>i = (Q!fe)fc>i, whence /„ = g n and r„ = s n = (j 2n -i/ {2n)) 2 /2. 

Step 4- We know that /„ is C°° and increasing on [0, r n ), that / n (0) = and that 
]im x ^.r n - fn(x) = 00 (due to (fT8|) and because j 2n > J2n-i)- Hence, there exists a 
unique g„ € [0,r„) such that f n (<ln) = 1 and we have < f n {q n ) < 00. Applying 
din]) at x = q n (recall that /„ = g n ), we deduce that q n f' n {q n )/n+{2—l/n) = l+2g„, 
so that /£(g„) = {n/q n )[2q n -1 + 1/n]. 

Step 5. Put 7„ = 2ri\J2q n £ (0, 2ny/2r n ) = (0,j2 n -i)- Then f n (q n ) = 1 rewrites 
7«^2n(7«) = 2nJ2n-i(7„). We now prove that j 2 n-i ^ 7« < 3 2n - 

• First, using (fT4|) with k = 2n — 1, we get xJ 2n (x) = (2n — l)J 2n ^i(x) — 
x J2n~i( x )- Since 7n^2n(7n) = 2nJ 2n -i(7„), we find Jan-i(7n) +7n^2n-i(7n) = 0. 
Thus 7„ > i 2 „_i, because for < x < j / 2n _ 1 < hn-u J2n-i{x) and J 2n -i( x ) are 
positive. 

• We next show that 7„ < j' 2n . We already know that -f n < j 2n -i- Wc thus 
assume below that j' 2n < j 2n -i, because else, there is nothing to do. There holds 
(2n/ 7 n)[7^2„-i(7n)]/[7^2„(7n)] = 1, whence (2n/7 n )[log( 7l 2 » J 2 „( 7 „))]' = 1 by 
(|16| . Consequently, 

[log J 2 „(7n)]' = (7n/2n) - (2n/7„). 

We know by (JTTJ) that j 2 „ > 2n. Hence J 2 „ > and thus [log J 2n ]' > on [0, 2n\. 
Thus 7„ > 2n, because for x < 2n, we have [log J 2n (x)}' > and (x/2n) — (2n/x) < 
0. Hence (7„/2n) - (2n/-y n ) > 0. Thus 7„ < j 2 „, because for x £ (j 2 „, 72n-i), 
J 2 „(:e) < so that [log J 2n {x)]' < 0. 

We thus have checked that 7„ £ (j' 2n -i, j' 2n )- Using (fTTj) . wc deduce that 

(2n - 1) + 2- 1 / 3 |a' 1 |(2n - l) 1 / 3 < 7 „ < 2n + 2- 1 / 3 |a' 1 |(2n) 1 / 3 + C(2n)- 1 / 3 , 

so that 7„ = 2n + la'-Jn 1 / 3 + O(l). Recalling that q n = (7„/(2n)) 2 /2, we easily 
deduce that q n = (1 + |ai|n~ 2 / 3 )/2 + C^n" 1 ) as desired. □ 

We now have all the tools to give the 

Proof 0/ Theorem^- (i) . Fix n £ N. Due to Lemma|ni there is a unique q n > such 
that f n (qn) = 1 and f^(q n ) < 00. Applying Lemma [5l we deduce the existence and 
uniqueness of a solution to (E n ). □ 

Let us now give two weak forms of the equilibrium equation (E n ). 
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Lemma 7. For n G N, consider the unique solution {c k )k>i to (E n ). For any 
4>, "0 : N i— > R with at most polynomial growth, 

n 



(20) 
(21) 



m (c") 
2 

mo(c n ) 



E [^( fc - 

k,l>l 

E ^( fc - 



o - m - m c n k c? = - e w*o - w)i ( fc - ^ 



k>2 



i) m\ hew = i e t^( fc ) - ^ fc ( fc - i)^- 



fe>2 



fc,Z>l 

Proof. Let n G N be fixed. 

Step 7. We first check that there is u n > 1 such that J2 k>1 u„c k < oo. This 
allows us to justify the convergence of all the series in Steps 2 and 3 below. Wc 



know (see Lemma [5]) that c£ 



V/n(9n)- Hence for any u > 0, J2k>i 



k~n 



U C 



fn(qnU)/ (q n fn(q n )). Since the radius of convergence r n of the entire series /„ 
satisfies r n > q n by Lemma[BJ the result follows: choose u„ > 1 such that q n u n < r n . 

Step 2. We now prove J20j). Multiply ©-flU by <j>[k) and sum for k > 1. We get 



m (c™) 



fc-i 



2m (c n )y>(fc)c2 



fc>2 



i=l 



fe>l 



E^)(fe-lK-[m 2 (c")-l]0(l) 



fc>2 



0(fc) E^-' c"c" 



Furthermore, there holds 



.>! + <t>(j)) cf c], as well as [m 2 (c") - 1]0(1) 



But Y,k>2 Y\ n J Z^i=l H "-ft 

2m„(c«) = £ 

[m 2 (c") - mi(c n )]0(l) = £ fe>1 fc</>(l)(fc - l)c£. This ends the proof of (2U) 

5iep 5. To check IpT]). it suffices to apply ([20} to the function </>(k) = kijj(k) 
and to use that by symmetry, J2k i>A( k + l)ip(k + I) — kip(k) — lxjj(l)]c^c" = 

2Ek,i>i k ii>( k + l )-^( k Mcf. a 

Remark 8. For n G N, consider the solution (c£)fc>i to (E n ), write p% = 
c k I 'mo(c") and then d k = p k /k. Then one easily checks, starting from i2U\) , that 
the sequence (d2)k>i solves, for all (j) : N i— > M. with at most polynomial growth, 



\ E ^ ■ 



i)-m-m Md n k d? 



k.l>l 



k>l 



k(k-l) 
2n 



di 



and has total mass ^2k>i k( ^k = Hence {d k )k>i is an equilibrium state for a 
coagulation-fragmentation equation with multiplicative coagulation kernel (particles 
with masses k,l merge at rate kl) and where each particle with mass k breaks into 
k particles with mass 1 at rate fc(fc — l)/(2n). This explains the similarities between 
Theorem^ and the results found by Rdth-Toth 32 ; the (non size-biased) particles' 
mass distribution (p2)k>l the same limit, as n — > oo, as the size-biased particles' 
mass distribution of the model considered by Rdth-Toth, see [321 Eq. (15)]. Observe 
however that the fragmentation rate in |32j is rather (fc — l)/«, which thus differs 
from k{k — l)/(2n). What seems important is just that roughly, for n very large, 
only huge particles break down into atoms. 



We are ready to handle the 
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Proof of Theorem^- (ii) . First, we know from Lemmas [5] and [6] that 

( n)= J_ = 1 1 1_ 

m ° [C ' QufM n(2q n -l + l/n) n[K| n -a/3 + 0(l/n)] ~ KK/3 " 

Next we use {H with 0(fc) = -1. This gives 

(22) mo(c n ) = m 2 (c") + m (c")-2m l( c") 

n 

whence, since m\{c n ) = 1, 

n 2/3 

m 2 ( c «) = (n-l)m (c") + 2~— . 

I oil 

Theorem [2]- (ii) is established. □ 

We can now prove the convergence of the particles' mass distribution. 
Proof of Theorem^ Let us put, for fc > 1 , 

2 (2k - 2 

Pfc ~ 4 fc fc v fc - 1 

Using the Stirling formula, one immediately checks ([7]). Next, recall that the Cata- 
lan numbers, defined by d = (^^j f° r au * — ®> satisfy, see [31] 

£-1 



C Y^CC : .-• (*>1), 



j'=o 



5^C i a; < = (l->/r^ 45)/(2x), x € [0,1/4]. 



i>0 



Observing that pk = 2.4 k C\-i, we easily deduce that p& = gX!i=i PiPk-i f° r 
> 2, as well as Efc>iPfc =2^m C fc _!(l/4) fe = E,> C,(l/4)*]/2 = 1. 

We now check that (j6|) has at most one nonnegative solution. To this end, it 
suffices to show that (j6|) implies that p\ = 1/2 (because this will determine the 
value of P2 = (l/2)pf, oip3 = (l/2)[pip2 +P2P1] and so on). To this end, it suffices 

to write Pl = l-J2 k >2Pk = WV 2 )E fe > 2 Ei=i = 1 - ( 1 / 2 )(E fc >i Pk) 2 = 

1 - 1/2 = 1/2. 

It only remains to prove that lim„ Efc>i |p£ — Pk\ = 0. Since Efc>iPfc = 
J2k>iP k = 1 , it classically suffices to prove that lim„p£ = p^ for all k > 1. First 
of all, we observe from ((5]) and Theorem EJ-(ii) that 

= cj = m 2 (c»)-l ^ n 2 / 3 /K| 1 = 

Pl m (c n ) 2nm (c n ) ~ 2n/(|a 1 |n 1 /3) 2 Pl 

as n — > 00. Next, we work by induction on fc. Assume thus that for some fc > 2, 
lim n p™ = p; for i = 1, . . . , fc — 1. Then, using (j4j, we deduce that 

c n y k—1 ^ fc— 1 

Pl = ^) = i+(t -i)/n &^ "> 2 = Pk 

as n — » 00. This concludes the proof. □ 
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We now study the size-biased particles' mass distribution. We start with the 
computation of all the moments and deduce a convergence result. 

Lemma 9. For each n G N, consider the unique solution (cj?)fc>i to {E n ) and 
the probability measure \i n = ^fc>i ^ c k^n- 2 / 3 k on (0, oo). 

(i) For any i > 1, 

where the sequence [Mi)i>o is defined by Mo = I, Mi = l/\a'i\ and. for i > 1 ; 
M i+1 = 2\a[\J2(^M j M i _ j _ 1 . 

(ii) There is a probability measure /j on [0, oo) such that for any </> € C([0, oo)) 
with at most polynomial growth, 

/•oo /-oo 

^<j)(n~ 2 l 3 k)kcl= <t>(x)ii n {dx) n -T 4>{x)n{dx). 
k>i J ° Jo 

This probability measure satisfies, for all i > 0, Jp 00 x l fi(dx) = Mj. 

Proof. Applying (f2Tj) with <fi(k) = k l with i > 1, we easily get 

E (•) m i+ i(c")m i _ j (c") - i K +2 (c") - m i+1 (c") - m a (c n ) + 1) , 



m (c") j-= 



so that 



2n * — 1 /'\ 

m 4+2 (c") = m i+ i(c n ) + m 2 (c") - 1 + ^— V ( * ) mj+i^Jmi-^c"). 

From this and the fact that we already know that 7712(0") ~ Afi?i 2 / 3 and that 
mo(c") ~ n — 1 ' /3 /|a' 1 1, one can easily check point (i) by induction. 

Next, it holds that J °° x l n n (dx) = ?i~ 2l / 3 ?77.; + i(c n ) for all i > 1. Consequently, 
we know from point (i) that lim„_ i . 00 x l /i n (dx) = Mj for any i > 1 . In partic- 
ular, sup n J °° x/j, n (dx) < 00, so that the sequence (/i n )n>i is tight. Furthermore, 
any limit point /i is a probability measure on [0, 00) satisfying J °° x l n(dx) = Mi 
for alH > 1. 

But such a probability measure /j is unique: an immediate (and rough) induction 
using that \a\\ > 1 shows that Mi < (2\a' 1 \)' l i\ for all i > 1. This implies the 
fmiteness of an exponential moment for /1, so that /1 is characterized by its moments. 
As a conclusion, fi n goes weakly, as n — > 00, to /i. This shows point (ii) for all 
<j) : [0, 00) 1 ^ M continuous and bounded. The extension to continuous functions 
with at most polynomial growth easily follows from point (i). □ 



Let us now show that \i satisfies some equilibrium equation. 
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Lemma 10. Consider the probability measure /i on [0, oo) defined in Lemma\$h 
For all 4> G C 1 ([0, oo)) such that <fi and <j>' have at most polynomial growth, 



2\a[\ 



oo poo 



(x + y)-(f>(x) , 

l{y>0} + <P {X)l{y=0} 



V 

x[4>{x) - cj)(0)]n(dx). 



H(dx)n{dy) 



Proof. We consider <fi as in the statement and, for n G N, \i n = J2k>i ^ c k^n~ 2 / ;i k 
as in Lemma[9l Apply ([21) with ip(k) = (f>(n~" 2 / 3 k): 

2n- 2 / 3 x - 0(n" 2 / 3 (fc + /))-0(n- 2 / 3 fc) ^ n 
m (c") ^ n-VH kc ^ 

n 2/3 



n 

k>l 



J2 (cp{n' 2/:i k) - 0(n~ 2 / 3 )) (n~ 2 / 3 fc - n^kcl- 



Multiply this equality by n 1 / 3 . In terms of fx n , this can be written as 

2n"V3 p /■<» cjj jx + y)-^) 

H n (dx)^ n (dy) 



m {c n ) J J y 

[<f>(x) - 0(n- 2 / 3 )](x - n- 2 / 3 K(c&). 



Recall that mo(c n ) ~ n 1 ' 3 /|ai| as n — » oo, so that 



OO /"OO 



r (x + y) -4>{x) 
2\a x \ I I fi n {dx)fi n (dy) 



o jo y 
[^{x)-^n- 2 ^)]{x-n- 2 ^)ii n {dx) 

[4>{x) - 4>(0)}xn n (dx). 

To obtain the last equivalent, use that <fr is continuous, that sup n J °° x/i n (dx) < oo 
and that sup n \4>(x) — (f>(0)\n n (dx) < oo by Lemma EJ-(i) since <fi has at most 
polynomial growth. Define now the function y) — ([<fi(x + y) — <fi(x)]/y)l{ y> o} + 
<//(x)l{j, = o}> which is continuous and has at most polynomial growth on [0, oo) 2 . 
Since \x n does not give weight to 0, we have 



OO fOO 



(23) 2K| / / T(x,y)fi n (dx) l i n (dy) 

Jo Jo 



oo poo 



Y (x + y) -4>{x) 
=2|o x | / / n n (dx)n n (dy) 



[<p{x) - (/)(Q)]xfj, n (dx). 

Recall that for any tjj G C([0, oo)) with at most polynomial growth, there holds 
L°° ijj(x)iJ> n (dx) Tl ^°° f°° ij)(x)(j,(dx) due to Lemma This implies that for all 
G C([0, oo) 2 ) with at most polynomial growth, J °° / °° ^(x, y) /i n (dx) n n (dy) 



16 



XAVIER BRESSAUD AND NICOLAS FOURNIER 



J °° J °° \& (x , y) fi{dx) n{dy) . Taking the limit as n — > oo in (|23|). we deduce that 

/>oo />oo />oo 

2\a[\ / T(x,y)n(dx)n(dy) = [<j>(x) - </>(0)]xn(dx) 
Jo Jo Jo 

as desired. □ 



We now try to determine a quantity resembling the Laplace transform of //. 
This is a usual trick for Smoluchowski's coagulation with constant kernel, see e.g. 
Deaconu-Tanre [T3] . 

Lemma 11. Consider the probability measure fi defined in Lemma\^ Then for 
all q > ; 

J x lailAi^ + o'J 

Proof. We put /3 = /i({0}), which we cannot exclude to be nonzero at the moment. 
We apply Lemma ITOl with <p(x) = (1 — e~ qx )/x, which is indeed in C 1 ([0, oo)) with 
(f>{0) = q and 0'(O) = ~q 2 /2. This yields 

2|oi|[A x (g) + A 2 (q) + A 3 (q) + A^q)) = B^q) - B 2 (q), 

where 

. , , A 00 A 00 + , 

Ai(q) = / / l {x>0 , y> o} V{dx)n{dy), 

Jo Jo V 



2{q)=p r m^m 1{y>oMdy) , 



o y 

oo 



A 

As{q) =/3 I (j)'(x)l {x>0} n(dx), 
Jo 

A 4 (q) =/?V (0), 
poo 

Bi(q) = I x<j>(x)fi(dx), 

oo 



II 



B 2 (q) =0(0) / Xfi(dx). 

Jo 

Recalling Lcmma[§]and that <fi(0) = q, we see that B 2 (q) = q/\a>i\- Next, 

/>oo 

Bi(g) = / (l-e-«*Mdx) = l-£'(q). 



One can check that <p'(0) = — q 2 /2, whence A±{q) = —j3 2 q 2 /2. A computation 
shows that [<fi(x) — <f>(Q)]/x + <j)'(x) = —q<fi(x), so that 

/•OO 

A 2 (q) + Ml) = ~H \ - e-^)\ {x> ^{dx) = -f3q[£(q) - f3q}. 

Jo 
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Finally, using a symmetry argument and then that (x + y)<p(x + y) — x4>(x) — y<p(y) 
-xy<t>(x)4>(y), 



oo />oo 



xMx + y) — x<j)(x) 

M(q) = / / l {x>0ty>0} fi(dx)fi{dy) 

Jo Jo X V 



oo />oo 



(x + y)4>{x + y)- X(j){x) - y<j>(y) 

Mx>o,y>o}V-{dx)n(dy) 

xy l a J 

1 ( r„ .-^ Kdx) 



2 



l-e\q)- q /\a[\, 



All this shows that 

2K| 
whence 

(24) ^(g) = l-g/K| + |oi|^(g). 

This equation, together with the initial condition £(0) = 0, has a unique maximal 
solution due to the Cauchy-Lipschitz Theorem. And one can check that this unique 
maximal solution is nothing but 

W \a[\M( q + a[)> 

which is dchncd for q £ (ai — a[, oo), because the Airy function does not vanish on 
(ai,oo). Indeed, it suffices to use that since Ai"(a;) = xAi(x), 

d_ ( Ai'{x) \ _ _ ( Ai'(xY " 

and that by definition, Ai^a'jJ =0. □ 

We now write down two formulae of Darling [TJ] and Louchard [35] that we 
found in the survey paper of Janson |26[ p 94]. Denote by [et)t£[o,i] the normalized 
Brownian excursion and define its area as *B CX = Jq 1 e t dt. Put, for y > 0, 

(25) lMy)=E[e-» B «]. 
There hold 

oo 

(26) ^ex(y) = v / 27 2 /^cxp(-2- 1 / 3 |a,|y 2 / 3 

(27) r ( l-e--)%S^ = 2V3^ i/ (°) Ai '( 21 ^) 



^2^3 W VAi(O) Ai(2V3g) 

This allows us to find a link between our probability measure /j, and tp ex . 

Lemma 12. Consider the probability measure (i defined in Lemma\^ Then 
ti(dx) _ yj cx (V2x^ 2 )e^\ 

~r ~ 20fk W 2 {x>0}dx - 
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Proof. Using Lemma ITT| we deduce that for all q > 0, 



= -^(K|) + %+K|) 

J_ / Ai^O) _ Ai'(g) \ 

ki Uko) Ai^;- 

Next, {27) implies that for all q > 0, 

1 /Ai'(O) Ai'(g)\ 
K\\M(0) M(q)J- 

We conclude by injectivity of the Laplace transform. □ 

We may finally give the 
Proof of Theorem^ Define c : (0, oo) i— > (0, oo) by 

(29) < X ) = ^r,, j 3 e /2 . 

By Lemma \\2\ the probability measure \i defined in Lemma |9] is nothing but 
fi(dx) = xc(x)dx. We thus have J Q xc(x)dx = 1. Recalling Lemma we know 
that lim„ ^ fc>1 (/>(n~ 2 / 3 fc)fcc^ = J °° (j>(x)xc{x)dx for any <f> £ C([0,oo)) with at 
most polynomial growth. Using (|26[) , we immediately deduce that 



(30) c(x) = |a / 1 |- 1 el a 'il x ^ 



-Kb 



It is clear from that c e C°°((0,oo)) and that c(x) ~ | a ' 1 |- 1 e ( :i 'il-l 0l l) a: : 
it suffices to use that < |ai| < |ci2 1 < that \dj\ (3nj/2) 2 / 3 and 



Lcbesguc's dominated convergence theorem. It is immediate from (|29l) that c(x) 
l/(2y / 7r|a' 1 |a; 3 / 2 ). Since we now know that n{dx) — c(x)dx does not give weight to 
zero, $5$i follows from Lemma [TU] Finally, © follows from Lemma [TT] when q > 0. 
It is easily extended to q > a\ — a\ using that c{x) ^2 \a[ | -1 ed a il-l a il) a \ □ 



4. Perfect simulation algorithms 

In this section, we provide some perfect simulation algorithms: we introduce a 
pruning procedure V n of trees (for each n > 1) that will allow us to interpret 

• the particles' mass distribution (p^)k>i as the law of the number of leaves of 
G n = VniG), where G is a binary critical Galton- Watson tree; 

• the size-biased particles' mass distribution (kc^)k>i as the law of the number 
of leaves of G n = Vn{G), where G is a size-biased binary critical Galton- Watson 
tree. 

We also show that G n obviously tends to G as n — >• oo, which gives a probabilistic 
interpretation of Theorem [3] 
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4.1. The pruning procedure. Consider a rooted discrete binary tree T, that is 
a set of vertices V(T) and of edges E(T), satisfying the usual properties of binary 
trees. The root is denoted by and we add a vertex * and an edge joining * to 0. 
We denote by L(T) the set of the leaves of T. We take the convention that * is not 
a leave, but may be a leave (iff V(T) is reduced to {*, 0}). We denote by N(T) 
the set of internal vertices (nodes) of T, that is N(T) = V(T) \ ({*} U L(T)). 

We recall that for any binary tree T, it holds that \L(T)\ = \N(T)\ + 1. 

Definition 13. Let T be a rooted discrete binary tree with at most one infinite 
branch and let n > 1. We define the (random) subtree V n (T) as follows. 

Step 1. Endow each edge e with the length n e , where (K e ) e ^E(T) * s an i-i-d. 
family of £xp(2) -distributed random variables. This induces a distance d on T . For 
each t > ; define T t — {x £ T : d(x, *) = t} and let to = inf{i > : T t = 0} < oo 
be the height of T . 

Step 2. For each internal node x £ N(T), consider the branch Bt(x) (endowed 
with the lengths introduced at Step 1 ) joining x to * and consider a Poisson point 
process it™ with rate 1/n on Bt{x) (conditionally on (ne)e£E(T)> these Poisson 
processes are taken mutually independent). All the marks of these Poisson process 
are activated. 

Step 3. We explore the tree from the top until we arrive at* ( that is, we consider 
T t for t decreasing from to to 0) with the following rule: each time we encounter 
an active mark of a Poisson point process (defined in Step 2), we remove all the 
subtree above the mark (and replace it by a leave) and we deactivate all the marks 
of the Poisson point processes corresponding to nodes in this subtree. See Figured 
for an illustration. 

Step 4- We call V n {T) the resulting tree, which is a subtree of T containing *. 

A mark of a Poisson process is said to be useful if has generated a pruning at 
some time in the procedure and useless otherwise. 

If T is infinite, then Step 3 looks ill-posed at first glance, since the top of the 
tree lies at infinity. 

Remark 14. Let T be a tree with one infinite branch. Then Definition Q3 
makes sense. Indeed, consider the first (starting from -k) edge e = (x\,X2) on 
this infinite branch, with xi child of x\, such that ir™ 2 has a mark in e. This 
happens for each edge independently with positive probability (not depending on e), 
so that such an edge a.s. exists. Then e will be cut (either by ir™ 2 or by a Poisson 
process corresponding to a node above X2 ) and this will make inactive all the Poisson 
processes corresponding to nodes above Xi . Thus everything will happen as if X2 
was a leave of T , neglecting all the marks due to Poisson processes corresponding 
to nodes above Xi- See Figured for an illustration. 

4.2. Particles' mass distribution. We can now give an interpretation in terms 
of trees of the particles' mass distribution. 

Proposition 15. Consider a binary critical Galton- Watson tree G (BCGWT 
in short), that is a Galton-Watson tree with offspring distribution (Sq + <52)/2. Fix 
n > 1 and let G n = V n (G), see Definition \13\ (all the random objects used by V n 
are taken conditionally on G). 
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(i) For all k > 1, Pr[|L(G)| = k] = pk, where pk was defined in Theorem^ 

(ii) For all k > 1, Pr[\L(G n )\ = k] = p^, where p 1 ^ was defined in Theorem^ 

Clearly, G n = V n {G) can be perfectly simulated. We thus have a perfect simu- 
lation algorithm for (pf)k>i- 

Heuristic proof. Any particle with mass k can be seen a cluster of k particles 
with mass 1. And it is natural to use a tree to represent the genealogy of this 
particle; if this particle has a mass k, then this tree will have k leaves. 




Figure 1 . The pruning procedure V n 

On Figure A, the tree T is drawn, its edges being endowed with i.i.d. £a;p(2)-distributcd 
random variables. The marks of the Poisson processes are represented as follows: here 7r™ 
has one mark (on the edge just under x), it™ has one mark (on the edge above on the left), 
7rJ has three marks (one on the edge just under z, one on the edge above on the right, 
one on the edge (*, 0)) and it" has one mark on the edge just under t. All the other Poisson 
point processes have no mark. 

Thus starting from the top of the tree, we first encounter the mark just under z. On Figure 
B, we have drawn the resulting tree: we have replaced the subtree above this mark by a 
leave and we have erased (deactivated) the marks sent by nodes of this subtree (here, the 
two other marks sent by z). 

Then we encounter the mark sent by t and the resulting tree is drawn on Figure C. Finally, 
we encounter the mark sent by x, which makes inactive the mark sent by y and the resulting 
tree is drawn on figure D. Since there are no marks any more, the tree of figure D is VniT). 
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To be more precise, we need first to handle some computations. Divide (|20j) by 
TO (c n ), use that J2k,i>i fiityPkP? = J2k>i HfyPk an d that, recalling ([22]) . 

i x> - d 2 0(ik = w ma(c " )+ 2 ( fir 2Trai(cm) = ^) = e 

?i * — ' nmn c" * — ' 

fc>2 uv ' k>i 

One gets, for all reasonable 0, 

2 + o - mWkp? + E^ 1 ) - mm 1 + ( fc - ^/^K = o. 

k,l>l fc>2 

Thus (p^)fe>i can be seen as the equilibrium of the mass of a particle with the 
following dynamics: (i) it merges with an independent similar particle at rate 1, 
(ii) its mass is reset to 1 at rate 1, (hi) its mass is reset to 1 at rate (k— l)/n, where 
k is its mass. 

Consider now such a particle at equilibrium. First neglect (iii) and follow the 
history of the particle backward in time: it has merged with a similar particle at 
rate 1 and it has been reset to 1 at rate 1. Thus this particle is subjected to events 
at rate 2 and each time an event occurs, it is a coalescence (node) with probability 
1/2 and a breakage (leave) with probability 1/2. This can be represented by a 
BCGWT, of which the edges have a length with law £xp{2). Now we take (iii) 
into account: we start from the past (i.e. from the top of the tree), we follow the 
branches of the tree and reset the particle to 1 at rate (k — l)/n, where k is the 
mass of the particle, i.e. the number of leaves of the subtree above the point under 
consideration. Using finally that the number of nodes of a binary tree is precisely 





FIGURE 2. Applying V n to a tree with one infinite branch 

The first edge e = (asi, X2) on the infinite branch such that 7r™ 2 has a mark on e is here the 
edge under t and X2 = t. Then independently of what can happen above f, this edge will be 
cut and all the marks sent by nodes above t will be erased. Indeed, there are two possibilities, 
calling M the mark sent by t. 

Case 1: This mark M is useful and thus M is replaced by a leave and all the marks sent 
between * and M by nodes above t are erased. 

Case 2. This mark M is useless and then it is necessarily deactivated by a useful mark lying 
between M and t (sent by a node above t). We conclude as in case 1. 

In any case, it is not necessary to know what happens above t to conclude that, with this 
configuration, T n {T) will be the same as in Figure [T}D. 
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its number of leaves minus 1, we guess that V n (G) should indeed provide a perfect 
simulation algorithm for the genealogy of a particle at equilibrium. 

We now handle a rigorous proof. 

Proof. We start with (i), which is completely standard. Define qk = Pr[|L(G)| = fc] 
for fc > 1 and use the branching property: knowing that G is not reduced to the 
root, it can be written as G = {0} UG'U G" , for two independent copies G',G" 
of G. Hence conditionally on {\L(G)\ > 2}, L(G) = L{G') U L(G"), whence, for 
fc > 2, {\L{G)\ = k} = {\L(G)\ > 2} n U^{\L(G')\ = i, \L(G")\ =k-i} and thus 
9fc = I Si=i QiQk-i- Next, it obviously holds that q\ = 1/2. Since pi = qi and 
recalling ((HJ), it follows that = p^. for all fc > 1. 

We now check (ii). Using the branching property of G recalled above, it imme- 
diately follows from the pruning procedure V n that for fc > 2, 

fc-i 

{\L(G n )\ = k} = {\L(G)\ > 2} n \J{\L(G' n )\ = i, \L(G%)\ = k - i,A n }, 

4=1 

where A n is the event that no pruning occurs in the edge (*, 0) and where G' n and G' 7 [ 
are two independent copies of G n , independent of ({L(G) > 2}, (7r™|( <ri 0)) 2 , e jv(G))- 
First observe that conditionally on (G^,G"), the event A n occurs if for all x G 
N(G' n ) U N(G' 1 [) U {0}, ttJ has no mark in (*, 0). Recalling that the length of this 
edge is £a;p(2)-distributcd, that the Poisson processes 7r™ have rate 1/n and that 
here we have |iV(G^)U7V(G^)U{0}| = \N(G' n )\+\N(G'^\+l = \L(G' n )\+\L{G»)\-l 
Poisson processes, one easily deduces that 

Pr[A n | G' nl Gl\ 



(\L(G> n )\ + \L(G»)\-l)/n + 2- 
Put now r£ = Pr[|L(G„)| = fc]. From the previous study, we get, for fc > 2, 

fc-i 



(i + k-i - l)/n + 2' 



whence [2 + (fc — l)/n]rj? = X)i=i r ? r fc-i- Using the arguments and notation of the 
proof of Lemma [5]-Step 1 (with r = 1), we deduce that rj£ = a^(ri) k for all fc > 1. 
But we also know, since G n C G is a.s. finite, that J2k>i a k( r i ) h = J2k>i r k = 1- 
Recalling (fTTj) and Lemma |H1 we conclude that r" = q n , whence rj£ = a^(q n ) k 
for all fc > 1. By Lemma [SJ it also holds that p£ = c^/mo(c n ) — ctf,(q n ) h , which 
concludes the proof. □ 

4.3. Probabilistic interpretation of Theorem [3j This is not hard from Propo- 
sition rr51 Clearly, as n — > oo, the probability that G = G n tends to 1, implying 
that the law of \L(G n )\, i.e. (p%)k>i, tends to the law of |L(G)|, i.e. (pk)k>i- 

Indeed, we have G„ = G as soon as for all x € N(G), 7r™ has no mark on the 
branch Bq(x) joining * to x in G. Conditionally on (G, {K e ) e eN(G))i this occurs 

with probability U x eN(G) CX P (~S £ee£(G),ec.B G (x) K e) , whi ch a.s. tends to 1, 
since G is a.s. a finite tree. 
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4.4. Size-biased particles' mass distribution. We now interpret the size-biased 
particles' mass distribution in terms of a pruned size-biased Galton- Watson. 

Definition 16. Consider a family of i.i.d. binary critical Galton-Watson trees 
{G(i))i>\. We call size-biased binary critical Galton-Watson tree (SBBCGWT in 
short) the binary tree G with one infinite branch (called the backbone) as in Figure 
where for each i > 1, we plant G(i) on *i. 




= *1 *2 *3 *4 *5 *0 *7 *8 *9 *10 *11 



Figure 3. A size-biased binary critical Galton-Watson tree 



Traditionally, each G{i) is planted above or under the backbone at random, but 
this is absolutely useless for our purpose, since we never take into account any order 
on the vertices. See Lyons-Pemantle-Peres [30l Section 2] for some indications about 
the terminology size-biased. 

Proposition 17. Consider a SBBCGWT G as in Definition[TBi and fix n > 1. 
Let G n — V n (G), recall Definition \13\ (all the the random objects used by V n are 
taken conditionally on G). Then for all k > 1. Pr(|L(G„)| = k) = kc^, where 
( c fc)fe>i was defined in Theorem^ 

Since V n (G) can be perfectly simulated due to Remark [T4l this provides a perfect 
simulation algorithm for (kc^)k>i- It seems striking that V n (G) is a size-biased 
version of r P n (G), with G a size-biased version of G. However, the two notions of 
size-biased are quite different: P n (G) is a version of V n {G) biased by the number of 
leaves, while G is a version of G biased by the size of the population at generation 
n (with n — > oo). 

Heuristic proof. First rewrite (|2ip as 

£ m + 1) - m]kc%2p? + - ^{ k )]^-±kc n k = o. 

k,l>l k>2 

Thus (fccj!)fc>i can be seen as the equilibrium of the mass of a particle with the 
following dynamics: (i) it merges with an independent particle with law (p^)k>i at 
rate 2, (ii) its mass is reset to 1 at rate (k — l)/n, where k is its mass. 

Consider such a particle at equilibrium and first neglect (ii). Then obviously, we 
can represent its genealogy as a forest of (non size-biased) particles V n {G{i)), the 
length of the edges one the backbone (those between these particles) being £xp(2)- 
distributed. Then, take (ii) into account: start from the past and prune the edges 
on the backbone at rate (k — where k is the mass of the particle, i.e. the 

number of leaves of the subtree above the point under consideration, whence k — 1 
is the corresponding number of nodes. 
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What we really do is slightly different, since we prune the (non size-biased) 
particles and the backbone simultaneously but one can easily get convinced that 
this changes nothing. 

We now give some rigorous arguments. 

Proof. Consider the problem with unknown (t^)k>i (here (p^)k>i is given) 

fe-i 

(31) X>2 = 1 ' [(A-l)/n + 2]t2 = 2 5>?#-< (*>2). 

fc>i »=i 

Step 1. Define s£ = Pr(|£(G n )| = k) for k > 1. We show here that (s£) fe >i 
is a solution to (|3Tj). First recall from Remark IT4l that G„ is a.s. finite, whence 
T,k>i s k = L Wc introduce G(l) the BCGWT planted on *1 and G' the SB- 
BCGWT on the right of *1 (see Figurc[3]). Clearly for k > 2, we can write 

{\L(G n )\ = k} = U* =1 {|£(G n (l))| = i, \L(G' n )\ = k - i} n A n , 

where A n is the event that there is no pruning in the edge (*, *1), where G„(l) = 
■p„(G(l)) and G' n = V n (G'). Note that G(l) and G' are independent and pruned 
independently so that G„(l) and G' n are independent. We know by Proposition 
ITSIthat |L(G„(1))| is (pj!)fc>i-distributed. Furthermore, it is clear that |L(G^)| has 
the same law as \L(G n )\. Finally, exactly as in the proof of Proposition [15] one 
may check that 

PrL4„ | G n [l),G' n ] = ^ . 

(|i(G„(l))| + |L(G;)|-l)/n + 2 

As a conclusion, there holds 

fe-i 

s k — Z^ Pi k ~ 



i=l 



(i + k-i- l)/n + 2' 



whence [(k — l)/n + 2]s£ = 2J2i=i Pi s k-i as desired. 

Step 2. Next, we show that (kc^)k>i also solves ([31]). Recall that ^fe>i ^ c fc = 1 
due to ©. For k > 2, using (gj), 

[(* - D/n + W = E ^^'^f^ = 2 g (fc ~yf% 

™o(c") ^ m (c") 

so that [(fc - l)/n + 2)kc n k = 2EtiP?(* - 

5£ep 5. To conclude the proof, it remains to prove that ([3T]) has at most one 
solution. For each x > 0, there is obviously a unique solution {u^(x))k>i to 

fe-i 

u?(x)=x, [(fc-l)/u + 2]^) = 2^p>LW (*>2). 

i=i 

One immediately checks recursively that for each k > 1, x i-> u^(x) is increasing. 
Hence, there is at most one value x n > such that J2k>i u k ( x n) = 1- But an y 
solution {t'k)k>i to (|3Tj) has to satisfy (t^)k>i = ( u fe(*i))fe>i an( l tnus must be 
equal to (u%(x n )) k >i. □ 
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5. A SCALING LIMIT FOR THE SIZE-BIASED TYPICAL PARTICLE 

The aim of this section is to understand why the Brownian excursion arises in 
Theorem 21 We will build a random real tree H and a pruning procedure Voo such 
that the measure of leaves of = Voo(H) follows the law xc(x)dx, where c is the 
profile introduced in Theorem [4] We will also explain heuristically why Hqo should 
be the limit of n -1//3 G n , where G n is the pruned SBBCGWT as in Proposition \17\ 
thus providing an interpretation of Theorem^ A rigorous proof of this scaling limit 
could probably be handled, but this would be quite tedious and not very interesting. 
For precisions about continuum random trees, see Aldous [2|4j[5], Duquesne-Le Gall 
[19] , Le Gall (2TJ [28] . Let us mention that some pruning procedures of continuum 
random trees have been introduced and studied by Abraham-Delmas-Voisin [l] , but 
it is quite different from ours. 

Wc proceed as follows: in Subsections 15.11 [5~2l and [5731 we recall basic facts about 
real trees, Aldous's continuum random tree and Aldous's self-similar continuum 
random tree. The pruning procedure is explained in Subsection [574] and wc indicate 
heuristically why the rescaled pruned SBBCGWT should converge to the pruned 
self-similar continuum random tree in Subsection 15.51 Subsection 15.61 is devoted to 
the proof that the number of leaves of the pruned self-similar continuum random 
tree is xc(a;)da;-distributed, where c was defined in Theorem [4] In Subsection 15.7] 
we briefly indicate how the pruned self-similar continuum random tree can be built 
in terms of contour functions. Finally, we discuss some remarkable property of a 
related diffusion process in Subsection 15.81 

5.1. Real trees. First recall the following definition, sec e.g. Le Gall [28] . 

Definition 18. 1. A measured rooted real tree T = (T, d, /i) is a metric space 
(T,d), with a distinguished element p called root, endowed with a nonnegative mea- 
sure \i, satisfying the following properties. 

(i) For all x,y £ T, there is a unique isometry (f>x,y '■ [0, d(x, y)] i— > T such that 
(f>x, y {0) = x and <j) x , y (d(x,y)) = y. 

(ii) For all x, y £ T, all continuous infective map h : [0, 1] i— > T with h(0) = x 
and h(l) = y, there holds h([0, 1]) = (j) Xty ([0, d(x, y)]). 

2. We denote by T the set of all measured rooted real trees. 

3. For T = (T,d,fi) € T, we denote by L(T) = {x £ T : T \ {x} is connected} 
the set of its leaves 

There is a classical way to build a compact real tree from a continuous nonneg- 
ative function, see Le Gall [28] . 

Lemma 19. Let a £ [0, oo) and g : [0, a] be continuous and satisfy 

g(0) = 0. For s,t £ [0,a], putd g (s,t) = g(s)+g(t)-2mm ue [ sAttSVt ] g(u), which is a 
pseudo-distance on [0,a]. We introduce the equivalence relation s ~ 9 t iff d g {s,t) — 
0. Consider the quotient space T g = [0, a]/ ^ 9 and the canonical projection tt 9 : 
[0, a] i y T g . The pseudo-distance d g induces a distance d g on T g . Define the root 
p = 7r 9 (0) and introduce the measure n g = A o tt^ 1 on T g , where X is the Lebesgue 
measure on [0,a]. Then T g = (T g ,d g , /i g ) is a measured rooted real tree in the sense 
of Definition \1S\ The function g is called the contour of Tg ■ 

This result is easily extended to build real trees with one infinite branch. 
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Lemma 20. Let g : R i— > R+ be continuous, satisfy g(0) = and g(—oo) = 
g(oo) = oo. For s,t € R wit/i s <t, we introduce 

( [s,t] tfO<s<tors<t< 0, 
KS ' ' J ~ \ ( — oo, s] U [t, oo) £/ s < < t. 

For s,t € R, pwi d g (s,t) = g(s) + .g(i) — 2 min ug ^ sAt sVt ^ g(it), which is a pseudo- 
distance on R. H^e introduce the equivalence relation s ^ 9 t iff d g (s,t) = 0. Con- 
sider the quotient space T g = R/ and £/ie canonical projection ir g : R i— > T s . 
TTie pseudo-distance d g induces a distance d g on T„. Define the root p = 7r g (0) and 
introduce the measure p g = A o 71^ 1 on T ff; where X is the Lebesgue measure on R. 
Then T g = (T g , d g , fi g ) is a measured rooted real tree in the sense of Definition \18i 
The function g is called the contour ofT g . 

5.2. Aldous's CRT. The famous continuum random tree (CRT) introduced by 
Aldous [31 SJ [S] is the real tree 72e = (T 2e , d 2e , M2e), where e = (e{) tg r 0j i] is a 
normalized Brownian excursion. The factor 2 is unimportant and replaced by 1 in 
many papers. Here we need the following one-parameter family of CRTs. 

Definition 21. Let e = (e t )te[o,i] be the normalized Brownian excursion. For 
u > 0, we introduce the Brownian excursion e" = (e") se [o. u ] with length u, defined 
by e" = ^/u.e s / u . We call CRT(u) the random rooted measured real tree T e ^ = 
(T e K , d e " , /i e " ) ( recall Lemma XW^l and we denote by Q u its law, which is a probability 
measure on T. 

The CRT is universal in that it is the scaling limit of any reasonable critical or 
sub-critical Galton- Watson tree conditioned on being large, as shown by Aldous [5]. 
Let us now recall some properties of the CRT we will use later. 

Lemma 22. Let u > and let T e ^ = (T e u , d e « , /i e « ) be a CRT(u) with root p. 

(i) p e u(T e u) = /i e „(L(T e «)) = u. 

(ii) J TgU d e u (p, x)p e u (dx) = J Q U e"ds = u 3/2 J Q e s ds. 

Proof. The fact that /z e «(T e « \ L(T e «)) = can be found in Aldous [5]. Now recall 
that by definition (see Lemma I19p , p e u = A o ir^} , where A is the Lebesgue mea- 
sure on [0,u]. Consequently, p e ™(T e u) = u. Furthermore, f Tu d e u(p,x)p e ™(dx) = 
Jq de" (p, 7r e « (s))ds. But p = 7r e «(0), whence, still by definition, d e « (p, 7r e u (s)) = 
d e . (0, s) = eft + e« - 2 min te [M ef = ejf . □ 

5.3. Aldous's SSCRT. The random real tree we will be interested in, because 
it arises as a scaling limit of the SBBCGWT, is, up to some scaling factor, the 
self-similar continuum random tree (SSCRT in short) of Aldous. The definition we 
give is similar to that of Aldous [3l Section 6], but the formalism is different. 

Definition 23. Consider [0, oo) endowed with its usual distance as a backbone. 

Let N(dx, du, dT) be a Poisson measure on [0, oo) x (0, oo) x T with intensity mea- 
sure 2dxn(du)Q u (dT) , where n(du) = du / X/2iru 3 . For each mark (x,u,T) of N, 
with T = (T, d, fi), plant T on the backbone at x (that is, connect the backbone and 
T by identifying the root ofT and x). At the end, we have a binary tree H , endowed 
with a distance 8 (induced by the distance on the backbone and the distances on the 
trees we have planted) and with a measure v (induced by the measures on the trees 
we have planted). Say that the root p is (of the backbone) . The random measured 
rooted real tree H = (H, 5, v) is called a SSCRT. It holds that v(H \ L(H)) = 0. 
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This last property is obviously inherited from Lemma [22]- (i). The following 
lemma corresponds to another construction similar to Aldous [4j Section 2.5]. 

Lemma 24. Consider a two-sided Bessel(Z)-process (Ct)tem, i-e. (Ct)t>o and 
{C—t)t>o ore two independent Bessel(3)-processes starting at 0, see Revuz-Yor [331 
Chapter VI] . The random measured rooted real tree 7c built from C as in Lemma 
\2(A is a SSCRT in the sense of Definition \2S\ 

All this is certainly contained in Aldous [3] [4] [5] with a different formalism. 

Proof. We divide the proof into several steps. 

Step 1. First, we use the excursion theory, see Revuz-Yor [33] Chapter XII], 
to build a Bessel(3)-process. For u > 0, denote by £ u the set of all continuous 
functions from [0, u] into [0, oo) and by R u the law of the Brownian excursion 
with length u, which is a probability measure on £ u . Set £ = U u> q£ u . Consider 
a Poisson measure M(dx,du,de) on [0,oo) x [0, oo) x £ with intensity measure 
dxn(du)R u (de). For x > 0, let t x = J Q J Q J £ uM(dx, du, de). The map x i— > t x is 
a.s. increasing on [0, oo). Define, for t > 0, its inverse L t = inf{x > : t x > t}. 
The map t H ► L t is a.s. continuous and nondccrcasing. Observe that for all x > 0, 
T x - = inf{£ > : L t = x} and t x = sup{i > : L t — x}. Consider X = {x > 
: r x - 7^ Tj;}. For x G I, we denote by (x, it(x), e(x)) the corresponding mark of 
M, with e(x) = (e(x, t)) 4e [ ,u(a;)] ■ By construction, u(x) = t x — t x - for all x G X. 
Define .4 = U re i(r x _, r^) and 6 = [0,oo) \ A. Put W t = for all t G B and 
Wt = e(x)e(x,i — t x _) for t G (ji-jTa), the sign e(x) being chosen with a fair coin 
toss for each x. It holds that IF is a Brownian motion and that L is its local time 
at and thus X = L + \ W\ is a Besscl(3)-process, see [33l Chapter VI, Corollary 
3.8]. It holds that L t = inf[ tj00 ) X s . 

Step 2. Build two Bessel(3)-processes X 1 and X 2 as in Step 1, using two inde- 
pendent Poisson measures M 1 and M 2 . Introduce the two-sided Besscl(3)-process 
C defined by C t = X\ for t > and C t = X 2 t for t < 0. Build the tree Tc as in 
Lemma 

Step 3. Introduce the Poisson measure N on [0, oo) x [0, oo) xT, image of M 1 +M 2 
by the map (x,u,e) i-> (x,u,T e )- The intensity measure of N is 2dxn(du)Q u {dT) , 
because the image measure of R u by the map e i— > T e is Q u , recall Definition 1211 
Build the tree H. = (H, S, v) as in Definition [23] with this Poisson measure N. 

Step 4- Finally, it is a tedious but straightforward exercise to show that 7c and 
% are a.s. isometric. Let us indicate that (i) the backbone [0,oo) of H corresponds 
{ttc(s) : (s > and L\ = X]) or (s < and L 2 _ s = X 2 S )} in T c (recall LemmaSD]) 
(ii) for each mark (x, u(x), e(x)) of M 1 the tree T e i x ) planted at x on the backbone 
of H corresponds to {itc(s),s G {t x _,t x )} in Tc (observe that C t = x + e(x,t — r x _) 
for t G {t x _,t x )), (iii) for each mark (x, u(x), e(x)) of M 2 the tree T e r x ) planted at 
x on the backbone of H corresponds to {7Tc(s),s G (— T 2 ,— t 2 _)}) in Tc (observe 
that Ct — x + e(x, — t 2 _ — t) for t G (— t 2 , —t 2 _)). We used the notation of Step 1 
with the additional superscripts 1,2. □ 

The following remark explains why the SSCRT should be useful to us. 

Remark 25. Let G be a SBBCGWT as in Definition UR For p > 0, endow 
the edges of G with a length p -1 / 2 . This induces a distance d p on G. Consider the 
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measure \x v = ^P 1 '^2 x( zg^x- Then we believe that (G,d p , p p ) tends to the SSC'RT 
(H,S,v) as p — > oo. A similar fact is mentioned in Aldous [21 Section 2.7] and 
Duquesne [181 Theorem 1.5] proves a very general result, including at least the fact 
that (G,d p ) tends to (H,S). 

Let us handle some brief computations showing that the scales and constant 
factors are correct. We fix x > 0. 

• Let Ap C G be composed of all the vertices (including roots and terminal 
points) of the Galton- Watson trees planted on {*, *1, . . . , ★LVp^J }> see Fig ure 131 
Then fi p (AP) = 4|A£|/p and one easily checks that = 1 + 2^!^ \L(G(i))\. 
Recalling that |L(G(1))| ~ (j>k)k>i and using the same arguments as in the proof 
of Theorem El we see th at E[ e - q ^ G ^} = 1 - \/l - e~<i , whence E[e~^( A x)] = 
e -4g/ P ( 1 _ ^ 1 _ e -8 g /p)LVP»J , which tends to cxp(-2x^). 

• Let A x C H be composed of the part of the backbone [0,a;] and of all the 
trees planted on [0, x]. Then v(A x ) = J* J T uN(dy, du, dT), since the backbone 
has no mass and since for each mark (y,u,T) of N, with T = (T,d,fx), it holds 
that p(T) = u (because T is a CRT(u) and due to Lemma |2"2"1 (i)). Consequently, 
E[e- q ^ A ^\ = exp(~x J °°(l - e -«") x /2/(7ru 3 )du) = exp(-2x^/iq). We used that, 
performing an integration by parts and the substitution u = x 2 , 

<*> r c - - r - r - ^ 

5.4. The limit pruning procedure. We introduce the pruning procedure Voo. 

Definition 26. Let T = (T, d, /i) be a measured rooted real tree. We define 
Vao{T) = (Toojdoc^oo) as follows. 

Step 1. Consider a Poisson measure M(dz,dx) on {{z,x) : z G T,x £ Bt(z)} 
with intensity measure (l/l6)/i(dz)\B T ( z )(dx), where Bt{z) is the branch of T 
joining z to the root p, endowed with the Lebesgue measure Xb t (z)- For each t > 0, 
we define Tt = {x £ T : 8{x, p) = t} and we introduce the height to = inf {t > : 
Tt = 0} £ [0, oo] of T ■ We activate all the marks of M . 

Step 2. We explore the tree T from the top until we arrive the the root p ( that 
is, we consider T t , for t decreasing from t to 0), with the following rule: each time 
we encounter a point X issued from an active mark (Z,X) of M , we remove the 
tree above X (i.e. the whole descendance of X ), we replace it by a leave and we 
deactivate all the marks {Z',X') of M with Z' descendant of X (which means that 
X £ B T (Z')). 

Step 3. We end with a tree T x C T , that we endow with the distance d^ and 
the measure p^, restrictions of d and p to T^. We set Voo(T) = {T^, dao, Poo)- 

A mark of M is said to be useful if has generated a pruning at some time in the 
procedure and useless otherwise. 

The picture is very similar to Figure [T] The pruning procedure Voo is not clearly 
well-defined if the height of T is infinite. 

Proposition 27. Let H = {H,8,v) be a SSCRT. Then V X {H) (the Poisson 
measure M used by Voo being picked conditionally on H.) perfectly makes sense: 
there is a.s. no need to start the pruning procedure from the top of H. 
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Proof. Let thus U = (H, S 7 v) be a SSCRT. We call B its backbone. For z G H, we 
define p(z) = argmin{<5(x, z) : x G B}. In other words, p(z) = z if z G B and p(z) 
is the root of the planted tree to which z belongs otherwise. 

For each a G N, consider the subset [a, a + 1] of B. Let f2 a be the event on 
which M({(z, x) : z G H, x G Bh(z) 7 p(z) G [a, a + 1], p(a;) G [a, a + 1]}) = 2, and, 
denoting by (Zi,X±) and [Z^^X^ the corresponding marks, there holds Xi G £>, 
X2 G B and a < X\ < p(Z\) < X2 < p{Z%) < a + 1 (which makes sense since 
a, X\, p(Zi), Xi,p{Z-z), a + 1 all belong to B). 

First, we claim that on fl a , the mark (Zi,Xi) is a.s. useful, so that we can do 
as if the backbone was ending at Xi and ignore all the marks (Z,X) of M with 
p(Z) > Xi (i.e. ignore all the marks sent by trees planted at the right of Xi). 
Indeed, on Q a , separate two cases. 

• Either the mark {Z^X-i) is useful, so that the backbone will be pruned at 
X2, whence all the marks (Z',X') with Z' > X2 will be deactivated. Furthermore, 
we know from f2 a that there are no marks (Z',X') satisfying Xi < X' < p(Z\) < 
p(Z') < X^,. Thus (Z\,X\) will indeed be useful. 

• Or the mark at (Z2, X2) is useless, so that it has necessarily been deactivated by 
an useful mark (Z%,Xz) with X3 G (X2, p(Z2)) (necessarily sent by a tree planted 
at the right of a + 1), but then the backbone is pruned at X3 and we conclude as 
previously that the mark (Zi,X\) will be useful. 

Next, (r2 a )aeN is clearly an independent family of events and it holds that p := 
Pr(fi a ) does not depend on a. We will check that p > and this will conclude the 
proof thanks to the Borel-Cantelli Lemma. 

To prove that p = Pr(fio) > 0, we observe that fij n £1q ^ ^0 C where 

• is the event that M (Ai) = 1, where A\ := {(z, x) : z G H,x G Bjj(z), p(z) G 
[0, 1/2], p(x) G [0, 1/2]} and that the corresponding mark (Z\,X\) satisfies X\ G B: 

• fig is the event that M(A 2 ) = 1, where A 2 ■— {(z, x) : z G H, x G B H (z),p(z) G 
(1/2, l],p(x) G (1/2, 1]} and that the corresponding mark (Z 2 , X 2 ) satisfies X 2 G B; 

• Oq is the event that M (A3) = 0, where A3 := {(z, x) : z G H, x G Bjj(z), p(z) G 
(1/2, 1], x G [0,1/2]} (here [0,1/2] C B). 

Recall that conditionally on H, M is a Poisson measure on {(z,x) : z G H , x G 
Bh(z)} with intensity measure m(dz,dx) := (l/lQ)v(dz)\B H ( z ){dx). Hence, know- 
ing H, the events fij, 51q and 51q arc independent (because the sets Ai,A 2 ,A 3 are 
a.s. pairwise disjoint). Knowing H, M(A\) is Poisson(m(A 1 ))-distributed and it 
obviously holds that m{A{) G (0,oo) a.s. (because v, restricted to the set of the 
trees planted on the subset [0,1/2] of the backbone is a.s. a finite and positive 
measure and because a.s., swp Z £H,p{z)£[o.i/2\ 

\B H (z)\ < 00). Hence, Pr(M(Ai) = 
1|H) = m(Ai) exp(— m(Aij) > a.s. But knowing "H, that M(A{) = 1 and knowing 
Z\ (here we denote by (Z\, X\) the mark of M in Ai), Xi belongs to [0, p(Zi)) C B 
with positive probability (since X\ is uniformly distributed on the branch Bh(Z\)). 
All this proves that Pr(fiQ|"H) > a.s. Similar arguments show that Pr(fiQ|'H) > 
and Pr(nl\H) > a.s., whence Pr(n \H) > a.s. and thus p = Pr(fi ) > 0. □ 

5.5. Heuristic scaling limit. Consider the pruned SBBCGWT G n = V n (G) as 
in Proposition IT7l endow its edges with the length iiT 1 ^ , which induces a distance 
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S n on G„. Also introduce the measure v n = An 2 ^ 3 ^2 x€ q S x on G n . Our aim in 
this paragraph is to explain heuristically why we expect that (G„, 8 n , v n ) converges, 
as n — > oo, to the pruned SSCRT VooiTL) = (i?oo, ^oo, v<x>)- 

Consider p very large, endow each edge of G with the length p~ 1 / 2 (which induces 
a distance d p on G) and set fi p = 4p _1 J~) x q S x . Then we know from Remark 1251 
that (G,d p ,fj, p ) resembles the SSCRT T~L = (H,S,u). Next consider n very large 
(with n >> p 1 / 2 ) and the pruning procedures V n applied to G. Then for each 
internal node z £ N(G), the Poisson process 7r™ sends 

• mark with approximate probability 1 — \Bq(z)\p 1 ^ 2 /(2ra), where |Bg(«)| = 
d p (p,z) is the length of the branch Bq(z) with G endowed with d p . Indeed, in 
the original scales, we have a Poisson process with rate 1/n on the branch B^(z) 
which is composed of p 1 ^ 2 \Bq(z)\ edges with £xp(2)-distributed length and the 
expectation of £xp{2) is 1/2; 

• one mark with approximate probability \B^{z)\p 1 ^' 2 / '(2n) (and then this mark 
is uniformly distributed in Bq[z)); 

• two or more marks with very small probability. 

Hence we can say roughly that the internal nodes z £ N(G) that do send one 
mark at x € Bq(z) are chosen according to a Poisson measure on {(z,x) : z € 

^ 1/2 

N(G),x € Bq(z)} with intensity m p>n (dx,dz) := ^ r Y^ ye N(G) S y( dz ) X B s (z){dx). 
Since each leave is very close to an internal node and since there are (roughly) as 
much leaves as internal nodes in a binary tree, we expect that 

pl/2 p 3/2 

m p , n (dx,dz) ~ 5 y (dz)\ B ^ {z) (dx) = —-^ p {dz)X Bs { z ){dx). 

yea 

Thus choosing p = n / 3 , whence p 3 / 2 /(16n) = 1/16, (G, d n 2/3,/i n 2/s) resembles the 
SSCRT H = (H,d,u) and the pruning procedure V n of (G, d n 2/3, fi n a/s) resembles 
the pruning procedure Voo of H. Calling finally S„ the restriction of d p to G n 
(which consists in endowing each edge of G n with a length p 1 / 2 = n 1 / 3 ) and u„ the 
restriction of u p to G„ (which gives a weight 4/p = 4n~ 2 / 3 to each vertex of G n ), 
we deduce that (G„, S n , v n ) resembles T-Loo- 

5.6. The measure of leaves. The main result of this section is the following. 

Theorem 28. Consider a SSCRT H = (H, S, v) as in Lemma \23\ and build 
'Poo('H) = {H OQl <5oo, Voo) as in DeRnition \26\ (picking the Poisson measure M used 
by Voo conditionally on 7-L). Then the law of v 00 {H 00 ) is xc(x)dx, where c was 
defined in Theorem^ 

Adopt the notation of Paragraph [531 Write n~ 2 / 3 |L(G„)| ~ n~ 2/3 |G„|/2 = 
I/ ™(G„)/8, which should tend to ^ 00 (iJ 00 )/8. Since |L(G n )| ~ (kc%)k>i by Propo- 
sition [T7] and since Voo{Hoo) /8 ~ xc(x)dx by Theorem |2"51 all this is coherent with 
Theorem |4j 

We start with two technical lemmas, of which the proof lies at the end of the 
section. 
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Lemma 29. Recall \25\) . There holds 

(l 5^ -* 

Lemma 30. Let p : [0, oo) H- (0, 1] be continuous, decreasing, satisfy p(0) = 1 
and limoo p — 0. Assume that for all q > 0, 



(33) 1-P(q)= p(s)h(q,s)ds, 

Jo 

where, the function I being defined by {Pp, 

h(q,s) = \a[\[£(( q + s)/2) - l{s/2)] M y^^ + ■ 

Ai (q/2 + aj 

T/ien /or a// q e [0, oo), = £'(g/2). 

Next wc compute, in some sense, the mass of the pruned planted trees. 

Lemma 31. For each u > 0, denote by m u the law of the mass of the pruned 
CRT(u). In other words, pickT u = {T u ,d u , (jl u ) a CRT(u), prune it as in Definition 
HI set Voo{T u ) = (T Ut00 ,d Uj00 , fi u ,oo): m u is the law of (M u ,oc(T Ut oo) ■ Define the 
measure n,^ on (0,oo) by 

poo 

n-ooiA) = l m u (A)n(du), 



where n(du) = du/Vinu 3 as in Definition \23[ It holds that 
noo (du) = e^ u/3 ^ ex (u 3/2 /16)n(du) 
where tjj ex was defined in 



As will be clear from the proof, we are not able to compute the law m„ for a 
given value of u. We really need to choose u according to the measure n(du). 

Proof. We introduce a Poisson measure N(dx, du, dT) on [0,1] x (0, oo) x T with 
intensity measure dxn(du)Q u (dT), wc denote by {(x, u x , Tx)}xex the set of its 
marks, where Tx = {T x , d x , p x ) with root p x . Recall that for each fJ- x (T x ) = 

u x by Lemma [2"2"1 Conditionally on N, we prune T x for each x € X independently, 
using Too with a Poisson measure M x on {(z, y) : z £ T x , y € P>t x (z)} with intensity 
measure (l/\Q)n x {dz)\ BT ^ z ){dy) and denote by T x ,oo = {T Xt00 , d x>00 , /J, Xt00 ) the 
resulting tree. Then we consider 



/>1 />oo /> 

Z:=^u x = / / 

xeI J0 JO JT 

Zoo '• — ^ p X ,oc (T X ,oo ) ■ 



uN{dx, du, dT), 



Step 1. Here we prove that for any r € (— oo, oo), 

-rZ 



exp (-J (1 - e"™) n 00 (du)J e [0,oo]. 
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Conditionally on {(x, u x ), x G I}, the random variables p, x ,oo(T Xj00 ) (for x £l) are 
independent and m« x -distributed, so that 



E(e- rZ °°\{(x,u x )} xex ) = [] / 



1 pOO 



log / e rv m u (dv) 



■ exp 

io Jo 

Taking expectations, we get 

/ pi poo 

E( e - rZ ~) =cxp ( / dx / n(du 

\J0 Jo JT 



u m Ulc (dv) 



N(dx, du, dT) 



Qu{dT) (e log /o^ - 1^ 

(/•OO / poo > 

/ n(cftt) I / e~ rv m u (dv) — 1 



(pOO poo 
/ n(cftt) / m u (dv) (e~ r " — ! 

We have finished, since n 00 (<iz;) = J °° n(du)m u (dv) by definition. 
>S£ep J?. The aim of this step is to prove that 



Pr(Z = Z c 



)\{(x,u x )} xex ) =exp( f f f \og^ cx {u 3 / 2 /16)N{dx,du,dT) 

\J0 JO J7 



We have Z = Zoo iff Tx.oo = T x for all x € I. Conditionally on N, the events 
{T Xt oo = T x }xex are independent. Furthermore, {T x ,oo = Tx) = {M x has no mark}, 
which occurs with probability exp(— J T p x (dz)d x (z, p x )/16), because for each z G 
T x , the length of Bt x (z) is d x (z, p x ). But knowing u x , T x is a CFQ^w-c), whence by 
Lemma [22| J T d x (z, p x )p x (dz) can be written as (w x ) 3 ^ 2 J e s ds, for some normal- 
ized Brownian excursion. We deduce that Pr(7i,oo = T x \u x ) = ip ex (u x ^ 2 / 16). Fi- 
nally, by independence, Pr(Z = Zoo\{(x, u x )} xeX ) = Hxex ^(("z) 3 ^/ 16 ): which 
ends the step. 

Step 3. Next we prove that there is A > such that 

Pr (Z = Zoo | Zoo) = cxp (-XZoo) ■ 
This is the most delicate part of the proof. 

Step 3a: Branching property for the forest of CRTs. Here we recall some results 
found in Duquesne-Le Gall [20l Section 4.2] (in the more general case of Levy 
trees). First, we denote by T = U x ^xT x the forest of CRTs, endowed with the 
distance d and measure p induced by (d x , p x )x<£ii and by T the corresponding 
measured forest (with d(y,y') — oo if y G T x and y' G T x < for some x ^ x'). For 
y G T, define the height of y as h{y) = d(y,p x ), where p x is the root of the tree 
T x to which y belongs. For a > 0, we introduce T(a) = {y G T : h(y) — a}, by 
T([0, a]) — {y G T : h(y) < a}, by T([0,a]) the corresponding measured forest, and 
we define similarly T([a,oo)) and T([a, oo)). 

For a > 0, we define the point measure N a on T(a) x T as A/" := J2iej °~(xf,S?)> 
where {Sf)i & j a G T are the closures of the connected components of T \ T([0,a]) 
endowed with the inherited distances and measures and xf is the root of 5". 

We also introduce the tr-finite measure Q on T defined by Q(.) = J °° n(du)Q u (.). 
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Then, by [20], there a.s. exists, for each a > 0, a finite measure k a on T such 
that the following points hold: 

(i) fc° = and for all a > 0, supp k a C T(a), 

(ii) for all a > 0, conditionally on T([0, a]), the point measure Af a is Poisson with 
intensity k a (dx)Q(dS) , 

(iii) M (.) = / °°fc%)<fa. 

S'tep 3b: Branching property for the marked forest of pruned CRTs. First, we 
call T the set of marked measured real trees: its elements are of the form (S, A) 
for S = (5, ds, fJ's) G T and for A a finite subset of L(S). Define the pruned forest 
%a = (Tc»)<ioojMoo) with = U xe iT X)00 similarly as in Step 3a. When pruning 
T to build Too, we keep track of the final useful marks and denote by A^ C T(Too) 
the leaves of Too where a pruning has occurred (i.e. A x = LfT^) \ L{T)). We 
put Too = (Tid^oo), which is a forest of elements of T. We also define Too (a), 
Too([0,a]), Too([0, a]), Too([a,oo)), Too([a,oo)) for a > as in Step 3a. We denote 
by %o([0,a}) = (Too([0,a]), Aoo nToo([0,a])) and define similarly %o{[a, oo)). 

For a > 0, we define the measure Af a ,<x> '■= J2i£j a ^ \ x <* s a ) on ^o°( a ) x 
T, where {Sf 00 )i<zj a arc the (marked and measured) closures of the connected 
components of Too \ T»([0, a]) together with ({y}, 0, 0, {y}) (the tree {y} endowed 
with the distance 0, the measure and with one mark at y) for y G Aoo n T 00 (a) 
(which is a.s. of cardinality or 1), and where xf^ is the root of Sf^. 

We finally introduce the measure Qoo(-) = /q°° **-(du)Q u>00 (.) on T, where Q u ,oo 
is the law of the marked pruned CRT(w) (that is, consider a CRT(w), prune it, keep 
track of the useful marks, this gives a marked measured tree, with law Q Ui00 ). 

The aim of this step is to check that, denoting, for each a > 0, by the 
restriction of k a (which is a measure on T) to Too (which is a subset of T), the 
following points hold: 

(i) fc°3 = and for all a > 0, supp k^ C Too (a), 

(ii) for all a > 0, conditionally on Too([0, a]), the point measure Af a ,oo is Poisson 
with intensity k^(dx)e- a ^ s ^ s ^ 1(i Q OD (d§) (with S = (S,d s ,Hs,A s ) G f), 

(iii) ^oo(-) = jr k %o(-) da - 

First, points (i) and (iii) are directly inherited from from Step 3a. To prove 
point (ii), we fix the height a. We are going to decompose the pruning procedure 
according to the decomposition of the tree T (at height a) introduced at Step 3a. 
Consider the Poisson measure M = X^xgi ^ x 011 {( z 'U) '■ z ^ T,y E Bt(z)} (recall 
from the first paragraph of the proof that the measure M x is used to prune T x , so 
that for x ^ y, M x and M y have disjoint supports). 

Given T, we write M = M a + J2ieJ (Mi + MX) where M a is supported by 
{(z,y) : z e T([0, a]),y G Bt(z)} while, for all i G J a , M° is supported by 
{(z,y) : z G Sf, y G B T {z),h(y) > a} and M° is supported by {(z,y) : z G 
Sf, y G B T (z), h(y) < a} = {(z, y) : z G Sf, y G B T (xf)}. 

We start the pruning procedure of T by M (from the top of T). Clearly, we 
first prune each subtree 5? (for i in J a ) using only M i . Conditionally on T([0, a]), 
the marked pruned tree <Sf ^ we obtain is measurable with respect to the er-ficlds 
generated by Sf and M i . It now follows from the branching property (Step 3a, 
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point (ii)) and from the definition of Qoo that, conditionally on T([0, a]), the point 
measure X^eJ ^(xf s a ) ^ s P°i sson with intensity k a (dx)Cloo(dS). 

For all i € J a , we consider the event on which all the marks sent by the leaves 
of Sf "under a" have been deactivated by the pruning procedure "above a" , i.e. 

^ = {m({{z,y) : z G St^y e B r (z), %) < a}) = 0}. 

This event is measurable with respect to (the cr- field generated by) Sf, Mf and M°. 
For all i e J"„, we have P(fi? | T([0, a]), 5?^) = e ~ Q ^ (5? -~ )/16 (here we used that 
afj,(Sf oo)/16 is the integral of the intensity of M° ). It follows that, conditionally on 
T([0, a]), the point measure Af°,oo := Ste,7 ^ fi °^(x a 5° ) * s P°i sson with intensity 
k a (dx)e- a ^ s( - s ^ 16 Q 00 (d§). 

On the one hand, 7^o([0, a]) is measurable with respect to T([0,a]) and to the 
(useful) marks of M falling under a (i.e. the useful marks (z, y) with h(y) < a). On 
the other hand, for all i <G J a , the event guarantees us that M° has no (useful) 
mark under height a. Thus %o{[0,o]) is, conditionally on T([0,a]), independent on 
■A'a.oo- Consequently, the restriction := J2iej a ln f 1 {^^u PP k^}\ x a § ; ^ of 

Afaoo to the support of k^ is, conditionally on 7oo([0,a]), Poisson with intensity 

k^(dx)e- a f" s ^/ 16 Q 00 (dS). 

We further observe that if i € J a ,oo, then is fulfilled: if Sf^ was sending a 
mark below a, then either this mark would be useful and remove S° M or this mark 
would not be useful and would be removed, together with Sf^, by another useful 
mark. This shows that A/ 0)00 is nothing but Af^oo- Point (ii) follows. 

Step 3c: Independence. Here we prove that for any a > 0, Aoo n T oo ([0,a]) is 
independent of A 00 nT 00 ([a, oo)) conditionally on Too- This means that conditionally 
on the (not marked) pruned tree Too-, the marks above a are independent of the 
marks under a. 

We fix a > and we consider the point measure oo := Sigj ^(x? ,5? ,5?) 
on (a) x T x T. The same argument as in Step 3b shows that, conditionally on 
Too([0, a]), it is a Poisson measure with intensity k OQ (dx)Il(dS , dT)e~ a ^ T< ^ T \ where 
R(.) = J °° n(du)Ru(.), for R u the joint law of (Soo,S), where S is a CRT(u) and 
6>oo is the resulting pruned tree. 

Conditionally on %o([0, a]) and %o, ( x i,oo^ioo)i&Ja,o>, i s completely determined. 
For each i £ J a .oo, the conditional law of Sf knowing 7^o([0, a]) and To depends 
only on Sf^ and in particular does not depend on Aoo H Too([0,a]). Since finally 
Aoo n ^('[o, oo]) = U i&Jai00 L(Sf t00 ) \ L(Sf ), the claim follows. 

Step 3d: Homogeneity. The aim of this step is to check that for < a < b, 
conditionally on To and knowing that #(j4oo H (Too (a) U T 00 (b))) = 1, we have 
#(Aoo n Too(a)) = 1 with probability fc^oo) / '[k^jT^) + k^T^)]. 

We consider two independent copies Too and TL- We are going to apply the 
branching property at height a for the first one and at height b for the second one. 
We define the point measure 



A MEAN FIELD FOREST-FIRE MODEL 



35 



on (Too(o) U T^ip)) x T. Independence (between Too and 7^,) and the branching 
property (Step 3b, point (ii)) shows that, conditionally on T>o([0, a]) and T£o([0, b]), 
the measure M a ,b,oo is Poisson with intensity measure 

(k^(dx)e- a ^ s ^ + k^(dx)e- b ^)Q 00 (dS). 

Consider the trivial marked tree f = ({0}, 0, 0, {0}) G T consisting of one 
(marked) point (the root). Standard properties of Poisson measures show that 
knowing Af a ,b,oo ((Too (a) U X£,(&)) x {?}) = 1 and denoting by (X, r) the corre- 
sponding mark of N a ,b,oo i it holds that conditionally on %o([0,a]) and 7^>([0, 6]), 
X G (a) with probability fc^(T 00 )/[fc^(T 00 ) + fc'L ( Hcr c we use that for 

5 = t, we have e'^ 3 ^ = e -b^s(S) = But (^f) i s a mark of A^b,^ if and 
only if X G n T^a) or X € A' x n T^fc). We have shown that 

Pr[#(4x> n roo(o)) - i|#((A co n T^a)) u «, n !£,(&))) = 1, Too] 

^oo (To ) 

^Sd(Txj) + ^'00(^00 ) 

To conclude, we fix /i G (a, 6). As a consequence of Step 3c, given Too and 
{To = TL}, the marked tree %o([0,h]) U ^([/i, oo]) has the same distribution 
as Too = Too([0, h]) U To{[h, oo]). Hence, the conditional distribution given Too 
and {Too = T4} and {#((Aoo H T^a)) U fl T^(6))) = 1} of the unique 
element of (A^ nToo(a)) U (A'^ FiT^Q))) is the same as the conditional distribution 
given Too and {#((^00 n Too(a)) U ((Aoo n Tx^o))) = 1} of the unique element of 
^00 n (Too (a) U Too (&))■ The claim follows. 

S'iep 5e: Conclusion. Using Steps 3c and 3d, one could prove that conditionally 
on Too, the set Aoo is a Poisson point process on i(Too) with intensity proportional 
to ^00 ■ We will prove a slightly weaker result, which suffices to conclude. 

Let us denote by Roo(a) = #{^4oo H TqoQO, a])}, which counts the number of 
marks below the height a of Tx). This process is clearly cadlag, nondecrcasing, and 
has only jumps of height 1 (it is clear from the pruning procedure that there will a.s. 
never be two marks at the same height). Step 3c implies that conditionally on Too, 
for all < a < b, Roo{b) — Roo{a) is independent of <j{{Roo{h) : h G [0, a]}). Hence 
(i?oo(a)) a >o is a Poisson process with some intensity measure moo (da) (this random 
measure being (r(Too)-measurable). Next, Step 3d tells us that for any < a < b, 
Pr[Ai?oo(a) = l|Ai?oo(a) + Ai?oo(&) = l,T>o] = k^T^) /[k^T^) + k^T^)]. 
We classically deduce that m 00 {da) is of the form Xk% (T 00 )da, for some A > 0. 
Consequently, 

Pr(i?oo(00) = 0|Too) = cxp (~ X f ^{T^da^j = CXp(-A^oo(Too)) , 

where we finally used Step 3b-(iii). Observe now that by definition, Zoo = Moo(Too), 
that {Z = Zoo} = {T = Too} = {Aoo = 0} = {#00(00) = 0} and that Zoo is ct(Too)- 
measurable. We conclude that 

Pr(Z = Zoo|Zoo) = E[Pr(i?oo(oo) = 0|Too)|Zoo] - cxp(-AZoo) 

as desired. 

Step 4- Writing E{e~ rZ l{ Z =z ac }] & [0,oo] in two different ways using Steps 
2 and 3 and that Z is a({(x,u x ),x G T})-measurable, we deduce that for any 
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r € (—00,00), I(r) = J(r), where 



=E 



cxp 



I(r) 

J(r) :=E[exp(-rZ t 
Using that 



rZ 



log ip ex (u s/2 /16)N(dx, du, dT) 



I(r) =E 



cxp 



[ru - logV> cx (w 3/2 /W)]N(dx,du,dT) 



■ exp 



1 



-ru+logi/> x(" 3/2 /16) 



n(du) 



as well as Step 1, we deduce that for all r, 

n oc (du) = 

and thus that for all q, 



1 _ e -(r+\)u 



e -™Vexv 



,3/2 



/16)) n(du) 



(1 - e-3") n^du) = r (l - e- ( «- A )"V cx ( u 3 / 2 /16)) n(cfo). 



Choosing g = 0, we see that J °° (l — e A ' ll '0cx(u 3 ^ 2 /16)) n(dit) = 0. This can clearly 
be satisfied only for one value of A, and we infer from Lemma [29l that 8A = \a[\ 
(use the substitution u = 8z). We have shown that for any q > 0, 

(1 - e-i u ) n^du) = j°° (l - e - 9 "el a 'il"/Vcx(" 3/2 /16)) n(du). 

Differentiating this equality with respect to q and using the injcctivity of the Laplace 
transform, we conclude that un^^du) = uS a ^ u ^ 8 ip e ^(u 3 ^ 2 /16)n(du) as desired. □ 

The following observation is somewhat obvious in view of Lemma 1311 

Lemma 32. Let B = [0,oo) stand for the backbone. Let R(dx,du,dy) be a 
Poisson measure on B x (0, 00) x (0, 00) with intensity measure r(dx,du,dy) := 
2dxn 00 {du){u/ld>)e~ uy / 1& dy, where n 00 (dM) was defined and computed in Lemma 
WT\ We activate all the marks of R and explore B from 00 to with the following 
rule: each time we encounter a point T = X — Y , where (X, U, Y) is an active mark 
of R, we prune B atT and we deactivate all the marks (X' ,U' ,Y') of R such that 
X' > T . We end with a finite segment [0, r] c B. Then 



T := 



B JO 



^{x<r}R{dx, du, dy) 



has the same law as ^oo(-ffoo), for some pruned SSCRT V^o (Ji) = (H 00 ,S oc> ,i' oc> ). 

Proof. The main idea is just that we can prune first the planted trees and then 
the backbone. Consider a SSCRT H = (H, 5, v) as in Definition [231 built with a 
Poisson measure N(dx,du,dT) with intensity measure 2dxn(du)Q u (dT), of which 
we denote by {(x, u x , T x ), x € 1} the marks. Conditionally on N, prune T x inde- 
pendently for each x € I, and denote by T x ,oo = (T x ,oo, d x>00 , A*x,oo) the resulting 
tree. Consider the point measure O = 5Z x ei ^0,m* =o(t x 00)) on B x (0, 00). Due to 
Lemma I3T1 O is a Poisson measure with intensity 2dxn OD (du). 

Then we prune the backbone as follows. For each x € I, we consider a Poisson 
process ir x on [0,x) c B with rate /i£,(T£,)/16. We activate all these Poisson 
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processes. Then we explore B from oo to with the following rule: each time we 
encounter a mark T of an active Poisson process 7r x , we prune the backbone at T 
and we deactivate all the Poisson processes ir v with y > T. We end with a finite 
segment [0, r] C B. Then [0, r], on which are planted the (pruned) trees T x ,oo for 
x € [0, T]nX, is a clearly pruned SSCRT Too {T~L) = {Hoo, <$oo, ^oo) and its total mass 
is nothing but ^oo(-ffco) = J2 x &x Hx,oo(T Xl00 )l{ x<T } = f B ul {x<r} 0(dx, du). 

Finally, we observe that for each x G X, at most one mark of tt x can be useful, and 
this mark is of the form x—Y x , where Y x is exponentially distributed with parameter 
Ma;,oo(Ia; i0 o)/16 (actually, tt 3 ' has one (or more) mark if x— Y x > and no mark if x— 
Y x < 0). Finally, we find the pruning procedure of the statement, with the Poisson 
measure R = Y^xex ^O-Mx ob(t x ^j.k,) which has indeed the intensity r and for which 
^00(^00) = J2 xe x ^x,oo{T x ^)l {x<T } = f B r Q °° f °° ul{ x<T }R(dx,du,dy). □ 

We can finally give the 

Proof of Theorem [2R We consider the random variable T built in Lemma [32] We 
have to show that T/8 ~ xc(x)dx. We introduce p(q) := E[e _<?r / 16 ] for q > 0. 
Recall from Theorem [4] that 



By injectivity of the Laplace transform, we thus only have to check that for all q > 0, 
p(2q) = £'(q). To this end, we will use Lemma [BUI Since p is clearly continuous, 
decreasing, satisfies p(0) = 1 and p(oo) = (because T > a.s.), it suffices to prove 
that p satisfies . We thus fix q > for the whole proof. 

Step 1. Recall the construction handled in Lemma 1321 Denote by {(x,u x ,y x ) : 
the marks of R. We consider the event tl(q) on which all the marks (X, U, Y) 
of R such that X — Y E [—q, 0] are finally deactivated (at the end of the procedure). 

Step 2. We prove that Pr(il(q)) = p(q). By construction, all the marks (X, U, Y) 
of R such that X > r arc deactivated (because r prunes the backbone) and 
all the marks (X, U, Y) of R such that X < t are still active and are known 
to satisfy Y > X (because else, the backbone would have been pruned at X — 
Y < t). Thus Pr[n(?)|{ar,u a } ai gi J r] = II iein[o,T] Pr t^ > x + q\Vx > x,u x ] — 
rLexnlo r ) e -9 " 1 / 16 , because y x is £a;p(M a ;/16)-distributed (conditionally on u x ). 
We immediately conclude, since X^GZn[o r) Ux = ^' * na ^ P r [^(?)l{ a; ; u a:}xei, r ] = 
e -?r/i6 ) w hence Pr(fi(q)) = E[e-« r / 16 ] as desired. 

S'fep 5. Next, we introduce 



and we denote by (V, Uv,Yv) and (W, Uw,Yw) the corresponding marks. We also 
consider the event A on which all the marks (X, U, Y) of R such that X > V 
and X — Y e [0, V] are finally deactivated when applying the pruning procedure 
described in Lemma [32] only on (V,oo) (i.e. exploring B only from 00 to V+). It 
holds that 




V := inf {x > : i?({x} x (0, 00) x (ar, x + q j) > 0} , 
W := inf {x > : R({x} x (0, 00) x (0, x)) > 0} , 



(n(q)) c = {V < w}n A. 



Indeed 
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• On {V < W} n A, all the marks (X, U, Y) with X > V and X - Y G [0, V] are 
deactivated (thanks to A) and there are no mark (X, U, Y) with < X— Y < X < V 
(because W > V), so that the mark (V, Uy, Yy) will not be deactivated and it holds 
that V — Yy G [—q, 0] by construction, whence cannot be realized. 

• On A c , there is a useful mark (X,U,Y) such that 0<X-Y<V<X, 
which thus makes (V,Uy,Yy) inactive (as well as all the marks (X',U',Y') with 
X' > V). Since by definition of V, there is no active mark (X, U, Y) with X < V 
and I-Fe [—q, 0], we deduce that 0(g) is realized. 

• On {V > W}, we have to consider two situations. Either the mark (W, U\y, Y\y) 
is useful, whence it prunes the backbone at W — IV G (0, V) and thus makes 
(V, Uy,Yy) inactive (and it also deactivates all the marks (X' , U' , Y') with X' > V) 
and we conclude as previously that Q(q) is realized. Or the mark (W,U\y,Y\y) is 
useless, but then it has necessarily been deactivated by a useful mark (X,U,Y) 
such that < X — Y < W which also deactivates (V, Uy,Yy), since V > W. We 
conclude again as previously. 

Step 4- By an argument of invariance by translation, we see that Pr[j4|V] = p(V), 
recalling that Pr[Qo(s)] = p(o) by Step 2. Indeed, conditionally on V — v, A is the 
same event as £Iq(v) concerning the Poisson measure R v = J2 X £Xn{v oo) ^(x-v,u x ,y x )- 
Next, the events A and {V < W} are independent conditionally on V. Indeed, 
conditionally on V — v, the event A concerns only the marks {(x, u Xl y x ), x G 
If! (v, oo)} while {W < v} concerns only the marks {(x,u x ,y x ),x G If] (0,v)}. 
Using Steps 2 and 3, we deduce that 

1 - p(q) = E(p(V) Pt[V < W\V}). 

Step 5. Define if>(q) = J °°(l - e - 9 "/ 16 )n 00 ((iw). Here we show that 

(34) E{p(V)-Pv[V < W\V}) = 2 r p{v) [^{v + q)~^{v)]e- 2 ^ +V ^ [x)dx dv. 

Jo 

First observe that for v > 0, {V > v} = {R{I V ) = 0} and {W > w} = {R(J W ) = 0}, 
where I v = {(x,u,y) : x G [0, v], u > 0, y G [x,x + q]} and J w = {(x,u,y) : x G 
[0, w],u > 0, y G [0,x]}. Since I v n J w = for any values of v > 0, w > 0, we 
deduce that V and W are independent. Recall that the intensity measure of R is 
r(dx 1 du 1 dy) := 2dxn OC} (du) (u/ 16) e~ uy / 16 dy. First, Pt(V > v) = e~ r ^\ with 

f v r°° f x+q u 

r(I v ) =2 dx n^du) J dy—e^/ 16 

=2 dx n oo (di0[e-"*/ 16 -l + l-e-"^ +9 >/ 16 ] 
Jo Jo 

=2 f dx [ip{x + q)- ip(x)] . 
Jo 

Consequently, the density of the law of V is 

hy(v) = 2 [ll>(v + q)- 1p{v)] e~ 2 tt +V ^ dx+2 Io ^) dx . 

Next, we have Pr(TU > w) = e^ 1 '^^ with 

r(J w )=2j dx J n^du) J dy^e^^ 16 = 2 J dx^(x). 
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By independence, we deduce that Pr(VF > V\V) = exp ^—2 dxip(x)\ . Formula 

(fMf immediately follows. 

Step 6. We finally deduce that p satisfies (|33[) . First, we observe from Lemma 
ED and Theorem H that 

noo(dw) = Jg-c(u/8)du. 

This implies, recalling that £(q) = /"(l - e- qu )c(u)du, that = \a' 1 \£(q/2)/2. 
We thus find, using Steps 4 and 5, that 

l-p(?)= / p(^)KI [l{{v + q)/2)-i{v/2)]e~ W ^ S * " l{x/2)dx dv. 
Jo 

To conclude, it only remains to observe from ([9]) that 

- Mfr e (x/2)dx =cxd ( r + " M'jx/l + al) \ = A?((g + v)/2 + a[) 
y \J q Ai{x/2 + a[) J Ai 2 (g/2 + ai) 

This ends the proof. □ 

It remains to prove the technical lemmas. We start with the 
Proof of Lemma\2^ Let us call / the integral in the LHS. Write I = K + L, where 

(l - el">) i>U^)^, L = £ (l - fe <^>) 5^. 

First, using (f28|) (that can be extended to all values of q € (oi, oo) by analyticity) 
with q = a[ and that Ai' (a[) = 0, we get K = Ai'(0)/Ai(0). Next, we write 
L = hmg^.oo(£g —L*), where 



A" 



20Fz 3 / 2 ' 9 7 v /mv ; 20Fz 3 / 2 ' 

Using j32]), we get L\ = y/q and {28]) implies L 2 = Ai'(0)/Ai(0) - Ai'(g)/Ai(g). 
Consequently, / = Um ? ^. 00 ( A /g + Ai'(g)/Ai(g)). But Ai'(g)/Ai(g) = -^ + 0(1/?) 
asg-^ oo, see e.g. Janson [26[ Section 16], so that / = as desired. □ 

We finally handle the 

Proof of Lemma [3b\ We first prove that there is at most one function p satisfying the 
assumptions of the statement. To do so, we first observe that, since t is increasing 
by definition, see (j9|), h(q,s) > as soon as q > 0. Next, a simple computation 
shows that for all q £ (0, oo), 

POO 

n(q) := / \h(q,s)\ds 



[ A i 2 ((q + s)/2 + a [) 
Ai 2 (g/2 + a;) 

= zi r 

Ai 2 (q/2 + a[) Jo 

Consider now two functions pi,P2 satisfying the assumptions of the statement, set 
u = pi — pi. It holds that ||u||oo — u (<7o) for some go G (0, oo), because u(0) = 0, 
limoo u = and because u is continuous. Observe next that 



Ai'((g + s)/2 + oi)Ai((g + s)/2 + a'Jds = 1. 



Mloo = \u(qo)\ = 



u(s)h(qo, s)ds 



< «(?o)||w| 
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whence ||w||oo = as desired since k(<?o) < 1- 

We now prove that p(q) := £'(q/2) satisfies all the requirements. First, by 
Theorem 21 there holds p(q) = e~ qx / 2 xc(x)dx, where xc(x)dx is a probability 
measure on (0, oo), so that p is clearly decreasing, continuous and satisfies p(0) = 1 
and linioop = 0. To prove that (|33|) . it suffices to check that H(q) = I(q) + J(q), 
where H(q) := Ai 2 {q/2 + a[)(l - £'(q/2)), where 

I(q) :=|oi| / e'(s/2)£((q + s)/2)Ai 2 (( + s)/2 + a'Jds 

£'(s/2)Ai((q + s)/2 + oi)Ai'((g + s)/2 + ctjds, 
(here we used the explicit expression © of £) and where 

/>oo 

J(q) := - |oi| / £ / {s/2)£{s/2)Ai 2 {{q + s)/2 + a'Jds 
o 

£ 2 (s/2)Ai((q + s)/2 + a[)Ai'((q + s)/2 + a[)ds 



(here we used an integration by parts and that ^(0) = 0). Recalling (|24[) . we deduce 
that \a' 1 \e 2 (s/2)~£'(s/2) = s/(2|ai|)-l, whence I(q) + J{q) = K(q) + L(q), where 

K{q) :=- Ai((q + s)/2 + oi)Ai'((g + s)/2 + a[)ds = Ai 2 (g/2 + a[) 

J<3 

and where 

1 f°° 

L{q) :=TTT7T / sAi ((° + s)/2 + a[)Ai'((q + s)/2 + a[)ds 



2\a[ 

1 



2\a[\ 
1 

2K| 



Ai 2 {(q + s)/2 + a , 1 )ds 
Ai 2 (s/2 + a' 1 )ds. 



Gathering everything, we realize that it only remains to check that L(q) = M(q), 
where (here we use (|24|) and then the explicit expression (j9]) of £) 

M(q) :=-£'(q/2)Ai 2 (q/2 + a' 1 ) 

= [1 - q/(2\a[\) + |a' 1 K 2 ( Q /2)]Ai 2 (g/2 + a[) 

= - M 2 (q/2 + ai) + qAi 2 (q/2 + oi)/(2|oi|) - (Ai'(g/2 + ai)) 2 /KI- 

Since lirtioo M = lirrioo L = 0, it suffices to check that M'(q) = L'(q), i.e. that 
2|ai|M'(g) = Ai 2 (o/2 + a[). But, recalling that Ai"{x) = xAi(x), 

2|oi|M'(g) :=(-2|oi| + ?)Ai(?/2 + ai)Ai'(<?/2 + a' x ) + Ai 2 (g/2 + oi) 
-2Ai'(g/2 + a' 1 )Ai"(g/2 + a / 1 ) 
={q + 2oi)Ai(?/2 + oi)Ai'(g/2 + ai) + Ai 2 (g/2 + ai) 
- 2(g/2 + oi)Ai(g/2 + oi)Ai'(?/2 + oi), 

which equals Ai 2 (g/2 + ai) as desired. □ 
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5.7. The contour process point of view. We now give an alternative presenta- 
tion of the pruned SSCRT. It might seem simpler at first glance, but when trans- 
posing the whole section to the formalism below, the proofs become rather less 
transparent. 

Remark 33. Consider a two-sided Bessel(i)-process (Ct)tem as in Lemma \^4\ 
Consider a Poisson measure 0(dt,dy) on the half plane R x [0,oo) with intensity 
measure (l/16)dsdy. 

Start with (D t ) t€ ^ = (C t )teR and with II the set of all the marks (T,Y) of O 
that fall under the graph of (Ct)tes. (i- e - f or which Y < Ct )■ 

Travel from y — oo to y = ( that is, consider the horizontal lines with ordinate 
y for y going from oo to 0) with the following rule: each time we encounter a 
remaining mark (T,Y) G H, we update (D t ) t £(L,R) an d H as follows. 

(i) Consider the largest (for inclusion) subset [(s,i)] C R containing T (this 
notation was introduced in Lemma WO) such that [(s,t)] x {Y} is under the graph of 
(D t )tt£R (i.e. such that for all u G C u > Y). 

(ii) Remove the zone [(s, t)] x [0, oo) from the picture (in particular, this makes 
disappear all the marks (T',Y') G II with T' G [(s,t)]) and glue if necessary: 

(a) if [(s, t}] = [s,t] for some < s < t, remove everything in the zone and then 
translate the picture on the right of t by t — s to the left; 

(b) if [(s,t)] = [s,t] for some s < t < 0, remove everything in the zone and then 
translate the picture on the left of s by t — s to the right; 

(c) if [(s,t)] — (— oo, s] U [t, oo) for some s < < t, remove everything in the 
zone. 

We end with a compactly supported function (D t )t e [L,R\, with L < < R, satis- 
fying Dl = Dr. We finally build the continuous function (C^°) t( =[o.R-L] by setting 
C t °° = D t for t e [0,R] and Cf = D L _ R+t for t G [R,R — L]. It holds that 
C^° = C^_ L = 0. Build the tree l~c°° as in Lemma \lfA Then T~c°° *s a pruned 
SSCRT V 00(H) = (Hoo, Soo, /ioc). In particular, T := R — L equals v^Hoo), so 
that r/8 is xc(x)dx- distributed. See Figure^ for an illustration. 

Let us explain briefly why this builds a pruned SSCRT. Recall from Lemma l24l 
that the SSCRT H = (H, 5, v) can be built from (Ct)teK as in Lemma 1201 consider 
the pseudo-distance dc on R, the equivalence relation ~ c on R, the quotient space 
Tq = R/ ^ c ., the canonical projection ttq '■ R <— > Tc, the distance dc induced by 
dc on Tc, the measure [ic — A o 71^ 1 on Tc (where A is the Lebesgue measure on 
R), the root p = 7rc(0). Define J t = inf[ t oo ) C u if t > and J t = inf^oo t] C u if 
t < 0. Set B = {ir c (t) : t G R, C t = J t }. Finally, set TL = (H, S, v) with H = T c , 
5 = dc,v = l±c- 

We introduce A = {(t,y) : t <E R, < y < C t }. For t > and y £ (0,C t ), 
we introduce s = sup{u G (0,t] : C u = y} and set Jt{y) = ^c(s). For t < 
and y G (0, Ct), we consider s = inf{u G [t, 0) : C u = y} and set Jt(y) = ^c(s). 
Observe that for all (t,y) G A, 7t(y) G H is the point of the branch BH{^c{t)) 
with 5(p,ir c (s)) = y. 

Next, we build a Poisson measure M(dz, dx) on {(z, x) : z G H,x G Bh(z)} with 
intensity measure (l/16)v(dz)\g H ( z )(dx), where Bjj{z) is the branch joining p to 
z, endowed with the Lebesgue measure. Denote by {(Tj, li)}i>i the marks of O in 
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A. For each i > 1, let Z t = nc{Ti) and X t = j Ti (Yi). Then M := Ei>i ^iA) 
is a Poisson measure on {(z,x) : z € if, x € Bjj(z)} with the correct intensity 
measure. Indeed, M is the image measure of O (restricted to A) by the map 
f(t,x) = (ivcit), 7t(x)). Hence conditionally on W, M is a Poisson measure on 
{(z,x) : z £ H,x £ Bh(z)} with intensity measure (1/16)Aa o / , A^ being the 
Lebesgue measure on A. Using that v = A(_ 00i00 ) o n^ 1 by definition, one easily 
concludes that M is as desired. 






G 




D 







H 



Figure 4. Pruning the SSCRT: the contour point of view 

On Figure A, the two-sided Bcsscl(3)-process (Ct)tem. is drawn, as well as the marks of the 
Poisson measure 0(dt, dy) with intensity (l/16)dtdy (we have erased all the marks above the 
graph of C). We assume that for t outside the picture (both on the left and right), Ct never 
goes under Cj^ , where Pi = (Ti , Yi ) . This allows us to claim that the mark P2 is useful (this 
is obviously the case if Pi is useful, and this is also the case if Pi is useless, because then it 
has been erased, as well as all the marks outside the picture. Thus P2 is useful. On Figure B, 
we thus build the corresponding zone, which we erase in Figure C. In Figure D, we consider 
the highest remaining mark and the corresponding zone, which we erase in Figure E. The 
last mark is treated in Figures F and G. Since there are no more marks, we have finished. In 
Figure H, we exchange the left and right parts of remaining function. Build the tree with this 
contour function as in Lcmma ll9l What we get is a pruned SSCRT 'Poo('H) = (Hoo, Sea, Moo). 
In particular, V = iAjo(-Hoo), so that T/8 is xc(x)dx-distributed. 
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Consider a mark (T,Y) of O, and the associated mark (Z,X) = (ttc{T), 7t(V)) 
of M . Assume that at some time in the procedure, the mark produces a pruning. 

(a) From the trees point of view (see V<x>), we remove the whole descendance of 
X and deactivate (forever) the marks (Z',X') of M with Z' descendant of X. 

(b) From the contour processes point of view (see Remark [33]), we consider the 
largest subset [(s,t)} c K such that [(s,t)] x {Y} C A and remove the whole zone 
[(s,t)] x (0,oo) from the picture. 

One easily understands that (b) is exactly the same thing as (a), because the 
subtree of H consisting of all the descendance of X is nothing but {ttc(u) : u G 
[(s, f}]} and the marks (Z',X') of M with Z' descendant of X correspond to the 
marks (T',Y') of O such that T G [(s,t)]. 

At the end, we get a function (-Di)te[L,H] . Roughly, (-Dt)te[o,fll is the "upper 
contour" of our pruned SSCRT, and (-Dt)te[L,o] is the "lower contour" of our pruned 
SSCRT. We need to identify the backbones. Exchanging the left and right parts of 
i D t)te[L,R] to get the function (C^°) te[0tR _ L] and using (C2°)te[o,R-L] as a contour 
function is a way to do it. 

5.8. A noticeable diffusion process. We conclude the paper with a remark 
about a diffusion process related to pruned CRTs. This is related to Lemma I3T1 
Consider a Brownian motion (Bt)t>o- Then exactly as in Lemma 124} (|-Bt|)t>o can 
be seen as the contour function of a family of CRT(u), planted at x G [0, 00) with 
a Poisson measure M(dx, du, dT) with intensity measure dx(2TTu 3 )- 1 / 2 duQ u (dT) 
(here there is no backbone: the CRTs are disconnected). Prune each of these CRTs 
independently using Voo ■ From the contour point of view (as in Subsection 15. 7p , 
we can handle the following procedure: consider a Poisson measure 0(dt, dx) on 
K x [0,oo) with intensity measure dtdx/16. Treat independently each excursion of 
(\B t \) t >o: start from the top of the excursion, etc... We obtain a new continuous 
nonnegative path (Yt)t>o, which is the contour function of an infinite forest. 

We believe that we can write Y t = \X t \, where X t solves the equation X t = 
dWt — \a' 1 \£(X t /2)dt/2. The main idea is that when at height x, the process X t is 
as a Brownian motion, but the pruning procedure pushes it down, and we found 
the drift coefficient by handling formal computations. We do not want to prove 
such a result, because the paper is long enough, but we believe that the following 
proposition is quite noticeable. 

Proposition 34. Fix a > and consider the odd function (3 : K h-> R defined 



Let (-Bt)t>o be a standard Brownian motion, let (X t ) t > be the unique solution to 



Jo 

and let {t x ) x >q be the inverse of the local time (L t )t>o at of (X t )t>o- Then it 
holds that E[e~ ATx ] = exp(— xip(X)), where 



by 



(3(x) 



sign(j;)a 



Ai'(a|a;| + a' x ) 
Ai(a|ar| + a[) ' 




(35) 
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where "0 e x was defined in \25\). In particular if a = 2 1 / 3 then we have the noticeable 
property that the drift coefficient (3 coincides with the Laplace exponent ip of the 
inverse local time. 

What is surprising in this result is that the drift coefficient and the Laplace 
exponent of the inverse local time are two objects of completely different nature. 

Proof. First observe that f3 is C 1 on R and has at most linear growth (because 
— Ai'(a;)/Ai(x) ~ ^/x as — > +00, see e.g. Janson [26l Section 16]), so that the 
S.D.E. has a unique solution. Next, the second equality in (|35j) is immediately 
deduced from ([9]). 

Using the result of Pitman- Yor [31] Theorem 1] , the Laplace exponent ip of the 
inverse local time (t x ) x >q of (X t )t>o is given, for A > 0, by 



(36) ip(\) 



>'a(0) 

<Mo) ' 

where <f>\ is a decreasing nonnegative solution to 

(37) \d>{{x)=p{x)<i>' x {x) + \<i> x {x). 

(With the notation of Theorem 1], it holds that ?p ' + = ?/> 0,_ because j3 is 
odd: for (j)\( x ) a decreasing nonnegative solution to (|37[) . 4>\(— x) is an increasing 
nonnegative solution to ([3"?) ). 

Define a := a- 2 Ai(2Xa~ 2 +a' 1 )/Ai(a[) and b := a - 1 Ai(2Xa- 2 +a' 1 )/Ai(a' 1 ). We 
claim that the solution to (|37p satisfying 4>\ (0) = a and 0' A (O) = b is nonnegative 
and decreasing. This will prove that ip(X) = /3(2X/a 3 ) and end the proof. 

First, this solution is obviously decreasing (and thus nonnegative) on R_ , because 
0a(O) = a > and (f>' x (0) = b < 0, whence 4>' x (x) < for all x < 0. Indeed, assume 
by contradiction that xq = sup{x < : §'\(x) = 0} exists. Then <j)' x < on (xo,0), 
so that (f)\ > (f>\(0) > on (x , 0], whence = 2[3<f)' x + \<f)\ > on (xo, 0) (because 
f3 < on K_), so that (j)'\(xo) < ^a(0) < 0, whence a contradiction. 

We next observe that for x > 0, 4>\(x) = ip\(x), where <fi\(x) := a _2 Ai(ax + 
2Aa -2 + a' 1 )/Ai(ax + a\). Indeed, it obviously holds that (p\(0) = a and ^' A (0) = b 
(recall that Ai'^) = 0) and some tedious computations, using only that Ai"(y) = 
yAi(y), show that ip\ solves (|37|) on M. + : check successively that 

(i) = a^ 1 Ai(ax + 2\a~ 2 + a[)/Ai(ax + cQ + ip\(x)/3{x), 

(ii) <p'l(x) = 2\<p x (x) + 2cp' x (x)(3(x) + <p x (x)\fi'{x) + a 2 (ax + a[) - p 2 {x% 

(iii) (3'(x) = -a 2 (ax + a'J + (3 2 (x). 

It only remains to check that ip\ is decreasing on M + (it is clearly nonnegative). 
But we see from (i) above that <p A (a:) = ip\(x)[(3(x) — (3(x + 2Aa~ 3 )] < 0, since the 
map x i — y P{x) is increasing on R + (use e.g. the second equality in (|35|> ). □ 
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